THE    GCARS    SYSTEM 
FORTRAN    M    PROGRAMMERS     MANUAL 
PART    C  -  PROGRAMS    FOR     CONTOUR     MAPPING 


SEPTEMBER     1969  -   NUMBER    27 


BY 


A    K.    TURNER 


JHRP 

JOINT  HIGHWAY  RESEARCH  PROJECT 

PURDUE    UNIVERSITY     AND 

INDIANA    STATE     HIGHWAY     COMMISSION 


Digitized  by  the  Internet  Archive 

in  2011  with  funding  from 

LYRASIS  members  and  Sloan  Foundation;  Indiana  Department  of  Transportation 


http://www.archive.org/details/gcarssystemfortr6927turn 


Final  Report   (Part  III  -  C) 

THE  GCARS  SYSTEM  FORTRAN  IV  PROGRAMMERS   MANUAL 
PART  C  -  PROGRAMS   FOR  CONTOUR  MAPPING 


by 

A.   Keith  Turner 
Post-Doctoral  Research  Instructor 


Joint  Highway  Research  Project 

Project:  C-36-72A 

File:  1-6-1 


Prepared  as  Part  of  an  Investigation 

Conducted  by 

Joint  Highway  Research  Project 

Engineering  Experiment  Station 

Purdue  University 

in  cooperation  with  the 

Indiana  State  Highway  Commission 

and  the 

U.S.  Department  of  Transportation 

Federal  Highway  Administration 

Bureau  of  Public  Roads 


The  opinions,  findings  and  conclusions  expressed  in  this 

publication  are  those  of  the  authors  and  not  necessarily 

those  of  the  Bureau  of  Public  Roads. 


Not  Released  for  Publication 


Subject  to  Change 


Not  Reviewed  By 

Indiana  State  Highway  Commission 

or  the 

Bureau  of  Public  Roads 


Purdue  University 
Lafayette,  Indiana 
September  1969 


Final  Report  (Part  III-C) 

THE  GCARS  SYSTEM  FORTRAN  IV  PROGRAMMERS 
MANUAL  PART  C  -  PROGRAMS  FOR  CONTOUR  MAPPING 


To:   J.  F.  McLaughlin,  Director 

Joint  Highway  Research  Project 

From:  H.  L.  Michael,  Associate  Director 
Joint  Highway  Research  Project 


September  19&9 
File:   1-6-1 
Project:   C-36-72A 


Attached  is  Part  III-C  of  the  Final  Report  on  the  research  project 
"Evaluation  of  Numerical  Surface  Techniques  Applied  to  Highway  Location 
Analysis"  which  was  approved  on  July  1,  1967.  This  report  has  been 
authored  by  Dr.  A.  Keith  Turner  of  our  staff. 

This  report  describes  in  detail  and  documents  the  FORTRAN  computer 
programs  used  by  the  GCARS  System  to  prepare  contour  maps  of  the  three 
dimensional  models  of  the  individual  or  combined  highway  location 
factors.  Contour  maps  are  prepared  on  the  printer  as  part  of  the  normal 
job  output,  or  they  can  be  produced  as  pen  and  ink  drawings  on  the 
CALCOMP  plotter.  The  report  is  one  of  a  series  of  four  reports 
documenting  the  various  phases  of  the  GCARS  System. 

This  report  is  presented  to  the  Board  for  review  and  comment  and 
for  the  record. 


Respectfully  submitted, 

Harold  L.  Michael 
Associate  Director 


HLM:mz 


F.  L.  Ashbaucher 
W.  L.  Dolch 

W.  H.  Goetz 

W.  L.  Grecco 

M.  J.  Gutzwiller 

G.  K.  Hallock 


M.  E.  Harr 

R.  H.  Harrell 

J.  A.  Havers 

V .  E .  Harvey 

F.  B.  Mendenhall 

R.  D.  Miles 


C.  E.  Scholer 

M.  B.  Scott 

W.  T.  Spencer 

H.  R.  J.  Walsh 

K.  B.  Woods 

E.  J.  Yoder 


NOTICE 

These  programs  have  been  tested  and  are  "believed  to  be  free  of 
errors.      If  errors  are  located,   the  Joint  Highway  Research  Project 
would  appreciate  notification. 

Send  comments  to:- 

Joint  Highway  Research  Project 
ATTN:     Professor  R.   D.   Miles 
School  of  Civil  Engineering 
Purdue  University 
Lafayette,   Indiana     ^7907 


TABLE  OF  CONTENTS 

Page 

LIST  OF  FIGURES  iii 

ACKNOWLEDGMENTS  i v 

INTRODUCTION  C-   1 

a)  Machine-Contouring  Procedures  C-   1 

b)  Computer  Display  Devices  C-  3 

METHOD  OF  OPERATION  C-  U 

PRESENT  CAPACITIES,   OPTIONS,   RESTRICTION  C-  6 


a)  Size  of  Matrices  C-  6 

b)  Irregularly -Spaced  Data  C-  6 

c)  Input  Options  C-  7 

d)  Output  Options  C-  7 


METHOD  OF  COMPUTATION  C-   7 

a)  Program  MAPPER  C-   7 

b)  Subroutine  GRID  C-  9 

c)  Subroutine  PRTCONU  C-13 

i)     Contour  Options  C-l8 

ii)     Display  Options  C-l8 

iii)     Contouring  Irregular  Areas  C-l8 

d}     Subroutine  HYPSO  C-20 

e)     Program  CONTUR  C-2". 

PROGRAM  OPERATING  CHARACTERISTICS  C-23 

a}     Program  MAPPER  C-23 

b)     Program  CONTUR  C-25 

PROGRAM  LISTINGS  AND  DOCUMENTATION  C-25 

a)  The  Use  of  FORSTAT  Data  C-25 

b)  Listing  of  Program  MAPPER  C-3C 

c)  FORSTAT  Cross-reference  Listing  for  Program  MAPPER  C-36 

d)  Listing  of  Subroutine  GRID  C-39 
FORSTAT  Cross-reference  Listing  for  Subroutine  GRID  G-^3 
Listing  of  Subroutine  PRTCONU  C-U6 
FORSTAT  Cross-reference  Listing  for  Subroutine  GRID  C-53 
Listing  of  Subroutine  HYPSO  C-56 
FORSTAT  Cross-reference  Listing  for  Subroutine  HYPSO  C-58 
Listing  of  Subroutine  STATS  C-60 
FORSTAT  Cross-reference  Listing  for  Subroutine  STATS  C-6l 
Listing  of  Program  CONTUR  C-62 
FORSTAT  Cross-reference  Listing  for  Program  CONTUR  C-69 
Listing  of  Subroutine  SCAN  C-72 
FORSTAT  Cross-reference  Listing  for  Subroutine  SCAN           C-7'+ 


11 


Page 

p)  Listing  of  Subroutine  TRACE  C-77 

q)  FORSTAT  Cross -Reference  Listing  for  Subroutine  TRACE          C-79 

r)  Listing  of  Subroutine  CALC  C-82 

s)  FORSTAT  Cross-Reference  Listing  for  Subroutine  CALC           C-83 

t)  Listing  of  Subroutine  DRAFT  C-85 

u)  FORSTAT  Cross-Reference  Listing  for  Subroutine  DRAFT         C-88 

v)  Listing  of  Subroutine  FROJ  C-91 

w)  FORSTAT  Cross-Reference  Listing  for  Subroutine  PROJ            C-91 


REFERENCES 


c-9? 


LIST  OF  FIGURES 

Figure  Page 

1  GGARS  Contouring  Programs  C-  5 

2  Flowchart  for  Program  MAPPER  C-  8 

3  Flowchart  for  Subroutine  PRTCONU  C-lk 
k                     Parameters  Used  in  Developing  a  Printer  Map  C-15 

5  PRTCONU  Display  Options  --19 

6  Definition  of  Contour  Range  Specifications 

Utilized  by  PRTCOI&  C-21 

7  Map  of  Drainage  Basin  Produced  by  PRTCONU  C-22 

8  HYPSOMETRIC  Curve  Produced  by  HYPSO  C-P2 

9  Geometric  Relationships  Used  in  Hypsometric 

Analysis  C-2U 

10  Typical  Card  Deck  Setup  for  One  Analysis  Using 

Card  Data  0-26 

11  Typical  Card  Deck  Setup  for  One  Analysis  Using 

Data  Stored  on  Tape  3  C-27 

12  Sample  Output  -  Program  CONTUR  C-28 


LIST  OF  TABLES 
Table  PaBr- 

1         Map  Distortions  Produced  on  Computer  C-17 


ACKNOWLEDGMENTS 

The  author  is  grateful  to  the  Joint  Highway  Research  Project  of 
Purdue  University,  the  Indiana  State  Highway  Commission  and  the  Bureau 
of  Public  Roads  for  their  financial  support. 

The  author  wishes  to  express  his  sincere  appreciation  to  Mr.  John 
Steele,  Associate  Director,  Purdue  Computer  Science  Center,  and  his  staff, 
without  whose  aid  the  documentation  of  these  programs  could  not  have  been 
completed  in  the  time  available.  Thanks  are  also  due  to  Miss  M.  Wong, 
Mr.  R.  Herrick,  and  Mr.  T.  Lassen  for  drafting  the  illustrations  and 
assembling  and  checking  the  computer  print-outs  utilized  in  this  reoort. 

Professor  Robert  D.  Miles,  P.  E.,  directed  the  research  project  and 
reviewed  the  manuscript.  His  assistance  is  gratefully  acknowledged. 


C-l 


THE  GCARS  SYSTEM 

FORTRAN  IV  PROGRAMMERS  MANUAL 

PART  C  -  PROGRAMS   FOR  CONTOUR  MAPPING 

by 
A.  Keith  Turner 

INTRODUCTION 

The  GCARS  System  is  a  computer -aided  highway  location  system.  Man- 
machine  communications  or  "dialogs"  are  required  between  the  engineer  and 
the  computer  during  all  phases  of  the  procedure.  Machine-to-man  information 
exchanges  must  be  concise  and  readily  comprehensible  if  the  system  is 
to  operate  efficiently.   In  most  cases  such  information  exchanges  involve 
the  production  of  both  (a)  statistical  and/or  summary  data,  and  (b)  graphical 
displays. 

Since  the  GCARS  System  utilizes  three-dimensional  data,  procedures  for 
displaying  such  data  were  developed.  This  section  of  the  GCARS  System 
Programmers  Manual  describes  the  programs  and  subroutines  used  to  generate 
contour  maps  of  three-dimensional  data, 
a)  Machine-Contouring  Procedures 

Contouring  of  three-dimensional  data  may  be  performed  in  a  variety  of 
ways;  the  contouring  may  be  extremely  literal,  extremely  interpretative,  or 
a  compromise  between  the  two  (l8)  .   Interpretative  contouring  involves 
considerable  "art",  since  the  contour  lines  do  not  necessarily  honor  all 
the  data  points  equally,  but  reflect  general  knowledge  of  the  process  by 
the  person  doing  the  contouring. 

Most  computer -generated  maps,  including  those  used  by  GCARS,  give 
literal  contours  that  honor  the  supplied  data.  Some  attempts  have  been  made 
to  modify  this  literal  interpretation  to  favor  certain  predetermined 
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directions  (12).  Even  with  literal  contouring  procedures,  the  results 
obtained  are  dependent  on  the  method  used  to  interpolate  between  points  of 
known  value.  As  a  consequence,  computer-generated  maps  are  not  "unbiased", 
as  is  often  believed  by  the  unwary,  but  are  in  fact  strongly  and  consistently 
biased  according  to  the  interpretation  of  the  data. 

As  a  consequence,  the  computer -generated  contour  maps  developed  by  these 
routines  may  not  be  equally  suitable  for  all  applications.  They  are 
extremely  useful  as  "quick-look"  maps  to  give  the  user  a  general  idea  of 
his  data,  or  to  check  the  general  validity  of  the  data.   Furthermore, 
since  all  maps  are  made  consistently,  they  can  also  be  used  for  comparison 
purposes.  Contouring  several  data  sets  by  humans  inevitably  involves  come 
learning  processes  which  introduce  an  uncertain  and  subtle  bias  into  the 
results  rendering  comparison  of  maps  difficult.   In  contrast,  the  computer 
will  not  "learn"  from  one  problem  to  the  other,  but  will  handle  each  data 
set  uniformly. 

Contouring  involves  identifying  and  connecting  a  sequence  of  points 
having  a  certain  value.  It  therefore  involves  searching  the  data  for 
areas  having  these  elevations,  and  then  computing  the  location  of  the 
contour  line.  The  only  convenient  way  to  do  this  on  the  computer  is  to 
use  gridded  data  -  a  series  of  "elevation"  values  of  the  dependent  variable 
being  contoured  on  a  suitably  scaled  square  lattice. 

Computer  contouring  therefore  generally  involves  two  steps  : 

(1)  conversion  of  irregularly  spaced  data  observations  defined  by 
X,  Y,  and  Z  coordinates  to  a  suitable  matrix  of  Z-valuer. . 

(2)  interpolation  within  grid  squares  to  locate  the  contour  lines. 
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If  the  data  is  already  in  grid  form,  the  first  step  can  be  omitted. 
This  is  not  a  common  occurrence  however.  A  large  number  of  interpolation 
procedures  have  been  proposed  including: 

(1)  fitting  of  a  suitable  regression  model  to  the  data  and  contouring 
the  regression  surface. 

(2)  fitting  of  "local  polynonials " .  This  involves  the  computation  of 
a  series  of  small  regression  models  utilizing  nearby  data  points, 
and  then  using  these  regression  surfaces  to  estimate  the  value  at 
the  grid  points. 

(3)  use  of  weighted  moving -average  procedures.  These  compute  an 
average  of  a  certain  number  of  nearby  points  with  each  weighted 
according  to  its  distance  to  the  grid  point  being  evaluated. 

Weighted  moving -averages  are  generally  computationally  easier,  and 
therefore  faster,  than  the  other  methods.  Few  quantitative  comparisons  between 
the  various  methods  are  possible.  Various  workers  prefer  certain  techniques 
based  on  their  experiences.  The  weighted  moving-average  procedures  were 
used  extensively  by  GCARS.  Subroutine  GRID,  based  on  a  program  developed 
by  Professor  W.  R.  Tobler  at  the  University  of  Michigan,  uses  a  weighted 
moving-average  algorithm  to  determine  grid  values.  As  in  all  interpolation 
procedures  some  smoothing  of  the  data  results  from  its  use.  A  limited  amount 
of  testing  suggests  that  the  smoothing  is  normally  less  than  2%   of  the  data 
range,  rising  to  nearly  10$  of  the  range  only  with  extremely  erratic  data. 
b)  Computer  Display  Devices 

Graphical  displays  can  be  developed  by  computers  utilizing  the  standard 
character  printer  or  one  of  the  specially  designed  display  units  such  as 
incremental  plotters,  drafting  machines,  or  cathode  ray  tube  viewers. 
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The  use  of  the  "on-line"  character  printer  often  has  advantages, 
particularly  in  the  initial  stages  of  an  analysis.   First,  the  graphical 
output  is  produced  by  the  printer  as  part  of  the  regular  job  output  and  is 
available  concurrently  with  the  analysis.  Second  the  printer  is  economical 
to  operate  and  relatively  fast  (600  to  1200  lines  per  minute).  Third, 
because  the  printer  is  an  "on-line"  device  it  does  not  require  the  use  of 
magnetic  tapes,  which  are  needed  when  the  CALCOMP  plotter  is  used,  thus 
reducing  the  chances  of  human  error. 

Purdue  University  has  an  "off-line"  30-inch  drum-type  CALCOMP  plotter. 
Use  of  the  CALCOMP  plotter  is  justified  when  plots  of  superior  quality  are 
desired.  The  resolution  of  the  printer  is  limited  by  line  and  character 
spacing.  Since  these  are  unequal,  distortion,  may  occur;  although  this  can 
be  avoided  by  careful  planning  of  the  display  scale  (see  page  C-lU).   In 
contrast  the  plotter  has  no  distortion  and  a  resolution  governed  by  its 
step  size  -  0.01  inch.  Thus,  it  can  produce  ink  drawings  suitable  for 
inclusion  in  reports. 

METHOD  OF  OPERATION 

Figure  1  shows  the  interactions  between  the  various  routines  described 
in  this  part  of  the  Programmers  Manual.  Subroutine  FRTCON^  which  prepares 
contour  maps  on  the  printer,  and  Subroutine  GRID  which  interpolates  irregular 
to  gridded  data,  are  utilized  by  several  portions  of  the  GCARS  System  (13,  1^). 
They  may  also  be  used  to  prepare  maps  of  data  in  conjunction  with  Program 
MAPPER,  as  shown  in  Figure  1. 

The  CALCOMP  contouring  program,  Program  CONTUR  and  its  associated 
subroutines,  is  a  semi -autonomous  sequence  of  routines  which  accepts 
punched-cards  containing  the  gridded  values  of  the  surface  to  be  contoured. 
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FIGURE   I.    GCARS  CONTOURING  PROGRAMS. 
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Subroutine  STATS  is  a  utility  subroutine  which  computes  statistical 
measures  of  data  arrays.  The  routine  determines  the  maximum  and  minimum 
values  of  the  data,  and  computes  the  mean,  standard  deviation  and  variance. 
It  is  widely  used  in  the  GCARS  System. 

PRESENT  CAPACITIES,  OPTIONS,  RESTRICTIONS 

a)  Size  of  Matrices 

The  GCARS  System,  as  presently  developed,  has  been  standardized  to  v;ork 
with  50  x  50  matrices.  Neither  dimension  can  exceed  these  limits  without 
correcting  the  DIMENSION  statements  in  the  programs.   Changes  in  shape,  for 
instance  to  100  x  25,  can  be  accomplished  wi thout  increasing  the  overall 
storage  requirements . 

Subroutines  GRID,  PRTCON1*,  and  STATS  automatically  change  their  matrix 
dimensions  to  coincide  with  those  of  the  calling  program.  These  programs  have 
been  successfully  run  with  matrices  dimensioned  100  x  100.   However,  due 
to  limitations  in  the  width  of  standard  computer  paper,  PRTC0NU  cannot 
contour  matrices  with  more  than  6k   columns  in  a  single  pass. 

b)  Irregularly-Spaced  Data 

Subroutine  GRID  can  make  interpolations  onto  a  grid  for  a  maximum 
of  1000  irregularly  spaced  data-points  at  a  time.   If  more  than  1000  data- 
points  are  located  within  a  desired  area,  it  must  be  subdivided  into  a 
series  of  subareas  each  containing  no  more  than  1000  points,  and  the  data 
must  be  presorted  into  appropriate  data  subsets.   If  the  chosen  subareas  are 
overlapped,  the  gridded  values  for  each  subarea  can  be  combined  to  give  an 
overall  grid.   Program  MAPPER  is  capable  of  pre-sorting  the  data  for  any 
area  into  a  number  of  subareas  before  calling  GRID.   It  checks  to  see  if 
the  subarea  contains  less  than  1000  points. 
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c)  Input  Options 

Program  MAPPER  can  be  supplied  with  either  gridded  or  irregular  data 
on  either  cards  or  tape.  Program  CONTUR  requires  card  data. 

d)  Output  Options 

Subroutine  GRID  can  punch  out  gridded  values  or  store  them  on  tape. 
Many  display  and  contouring  options  are  available  in  both  PRTCONU  and  CGI7TUR. 

METHOD  OF  COMPUTATION 
a)  Program  MAPPER 

Program  MAPPER  has  been  developed  to  produce  contour  maps  on  the  printer 
for  either  gridded  or  irregularly  spaced  data.  The  maps  ^re  produced  by 
Subroutine  PRTCONU,  which  required  gridded  data.  Thus  if  the  data  is  not 
gridded,  MAPPER  first  calls  Subroutine  GRID  to  interpolate  the  irregularly- 
spaced  data  onto  a  suitable  grid. 

Figure  2  is  a  flowchart  of  Program  MAPPER.  The  program  reads  the  first 
control  card  to  determine  the  type  of  data  (gridded  or  irregular)  and  branches 
accordingly.   In  either  case  the  input  data  may  be  supplied  on  cards  or  on 
tape.  In  accordance  with  the  GGARS  standards,  irregular  data  is  supplied 
on  a  tape  labelled  TAPE2,  and  gridded  data  on  a  tape  labelled  TAFE3.   (An 
option  within  Subroutine  GRID  allows  the  user  to  request  the  interpolated 
grid  values  be  stored  on  a  TAPE3  user  tape). 

For  computational  efficiency  the  number  of  irregularly  spaced  observations 
being  interpolated  at  any  one  time  is  limited  to  one  thousand.   The  total 
running  time  of  Subroutine  GRID  depends  on  the  number  of  searches  and 
comparisons  it  has  to  make.  The  number  of  comparisons  is  the  product 
of  the  number  of  grid  points  being  evaluated  and  the  total  number  of  original 
data  points. 


FIGURE    2      FLOWCHART    FOR    PROGRAM    MAPPER. 
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Accordingly,  Program  MAPPER  allows  for  the  specification  of  several 
overlapping  subareas  and  "pre-sorts"  the  new  data  to  locate  those  points 
falling  within  each.  The  total  is  checked.  If  it  exceeds  1000  points  an 
error  message  is  printed  and  further  analysis  of  that  subarea  is  abandoned. 
The  resulting  output  will  consist  of  several  card  decks  and  printer  maps, 
which  can  be  combined  to  form  a  single  gridded  array  and  a  single  map  of  the 
whole  area  (Programs  for  this  step  are  left  for  the  user  to  develop).   Once 
a  single  array  has  been  produced  it  can  be  resubmitted  to  produce  overall 
maps  of  the  area. 

After  the  data  have  been  prepared,  the  two  branches  of  Program  MAPPER 
converge  and  Subroutine  PRTCONU  is  called  to  produce  contour  maps  of  the  data. 
Multiple  maps  of  the  data  can  be  produced.  These  may  be  different  or  identical. 
Tf  identical  maps  are  desired,  the  word  "SAME"  must  appear  on  the  appropriate 
control  card,  and  then  only  one  set  of  map  control  cards  are  required. 
Otherwise  the  number  of  map  control  card  sets  must  equal  the  number  of  maps 
requested.  Details  of  the  various  control  cards  are  given  on  page  C-30  as 
part  of  the  listing  of  Program  MAPPER. 
b)  Subroutine  GRID 

The  documentation  for  an  interpolation  algorithm  and  a  MAD  language 
program  were  obtained  from  Professor  W.  R.  Tobler,  Geography  Department, 
University  of  Michigan,  Ann  Arbor.  The  MAD  language  program  was  translated 
into  FORTRAN  IV  and  forms  the  basis  of  Subroutine  GRID  within  the  GCARS 
System. 

The  algorithm  is  used  whenever  irregularly  spaced  observation,  each 
defined  by  X,  Y,  and  Z-coordinates,  are  interpolated  to  produce  a  matrix 
of  X-coordinate  values,  corresponding  to  the  positions  of  a  series  of 
observations  on  a  regular  lattice. 
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When  using  Subroutine  GRID,  the  operator  may  define  the  limits  of  the 
area  over  which  interpolation  is  to  be  performed,  or  the  routine  can  search 
for  the  limits  of  the  area  defined  by  the  data  and  use  these  as  the  inter- 
polation boundaries.  The  actual  spacing  of  the  lattice  can  be  defined  by 
the  operator.  Otherwise  the  program  will  compute  a  "best"  grid  size  which 
approximately  equals  the  square  root  of  the  area  divided  by  the  number  of 
observations.  This  procedure  is  described  in  detail  below.  Since  the 
interpolation  boundaries,  the  grid  sizes  in  the  X  and  Y  directions,  end 
the  number  of  columns  and  rows  computed  or  defined  for  the  matrix  are  all 
mutually  inter-dependent,  some  readjustments  may  be  required. 

The  algorithm  involves  a  weighted  moving  average  procedure,  and  appears 
to  give  results  which  are  comparable  to  those  obtained  by  the  fitting  of 
local  polynomials,  but  is  computationally  much  simpler.  Following 
considerable  (empirical)  testing,  the  main  disadvantage  appears  to  be 
a  slight  dampening  (or  smoothing)  effect  on  the  data.  This  approaches  10$ 
of  the  data  range  for  extremely  erratic  data,  but  for  commonly  encountered 
data  distributions  is  much  less,  one  or  two  percent  smoothing  being  normal. 
Competing  schemes  are  described  in  the  references  (l,  k) . 

The  following  discussion  was  obtained  from  Tobler.   It  includes  a 

description  of  the  procedure  for  defining  an  optimum  grid  size. 

An  observation  is  recorded  as  X. ,  Y. ,  Z. ,  i  =  1,  ...  ,  n,  where 

i'  1   1'  '   ' 

X,  Y  are  locational  coordinates  and  Z  is  a  value  at  X,  Y.  Lattice  points 
are  labelled  jk  in  standard  matrix  notation  (j  =  row  subscript,  and 
k  =  column  subscript).  The  problem  is  to  estimate  Z1   . 
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This  is  accomplished  in  the  f ollowing  steps : 

(1)  Calculate  A,  the  area  of  the  rectangle  bounded  by  the  maximum 
and  minimum  of  the  coordinate  values  X  and  Y. 

(2)  Let  the  grid  size,  G,  equal  the  square  root  of  the  area 
divided  by  the  number  of  observations, 

G  =  {k/nf 

(3)  The  number  of  rows  in  the  matrix  is  computed  as  the  (truncated) 
integer  part  of 

R  ~   (A  Y/G)  +  1.5  where  A  Y  =  range  in  Y 
and  the  number  of  columns  as  the  integer  part  of 

C  =  (A  X/G)  +1.5  where  A  X  =  rangp  in  X 
This  yields  a  matrix  which  has  approximately  N  entries. 
(h)     In  general,  because  of  truncation, 

(C  -  l)  *  G  /  A  X 

(R  -  1)  *  G  /  A  Y 
Split  the  difference,  letting 

A'  X  =  (AX  -  (C  -  1)  *   G)/2 

A'  Y  =  (AY  -  (R  -  1)  *  G)/2 

and  assign  X  and  Y  coordinates  to  positions  in  the  matrix. 

For  example  the  j   row  and  the  k   column  will  have  the  coordinates 

X,  =  X  ,   +  A'  X  +  (k  -  l)  *  G 
k    min         v     ' 

Y .  ■  Y    -  A'  Y  -  (j  -  1)  *  G 

j    max         vu    ' 

(5)  For  each  position  in  the  matrix  calculate  the  square  of  the 
distance  to  each  observation  using 
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and  rank  the  points  in  order  of  increasing  distance  from  the 

jk   lattice  point  as  they  are  computed,  saving  the  results 

2 
for  only  the  nearest  six  points,  except  if  D.   <  e  in  which 

case  set  Z'   =  Z.  and  proceed  to  the  next  position  in  the  matrix 

JK     1 

since  the  assignment  has  been  completed,  (e  has  been  arbitrarily 
set  =  G/25.  This  means  that  if  any  data  point  is  as  close  or 
closer  than  l/25  of  the  grid  size  to  a  lattice  point,  the 
lattice  point  is  assumed  to  be  adequately  defined). 
(6)  Assign  a  value  to  the  jk   position  by  computing  the  weighted 
mean  as  follows : 

(a)  compute  the  weighted  mean  of  the  nearest  six  points 

as  the  sum  of  the  elevations  divided  by  the  square  roots 

t 
of  their  distances  to  Z   ,  divided  by  the  sum  of  one 

over  the  square  roots  of  these  distances. 

(b)  average  this  weighted  mean  with  the  value  of  the  nearest 
of  the  six  points.   The  algorithm  may  be  written  as 


7          + 

6 

Z.    z./i/d. 
iti     l/  V    1 

n 

i  1/  xfX 

Z'-v 


where  Z'  .,  is  the  desired  interpolated  value 

Z    is  the  nearest  data  value 
n 

Z.   are  the  six  nearest  data  values 

1 


D.   are  the  corresponding  distances  from  Z.  to  Z'  .,  . 
1  f  &  1      jk 
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c)  Subroutine  ERTCOM^ 

This  routine  is  also  derived  from  a  MA.D  language  program  obtained  from 
the  Department  of  Geography,  University  of  Michigan,  Ann  Arbor.  The 
program  accepts  a  matrix  of  elevation  values.  These  are  assumed  to  form  a 
square  lattice  (i.e.  the  spacing  in  the  X-direction  between  columns  equals 
the  spacing  in  the  Y-direction  between  rows ) . 

Figure  3  is  a  flowchart  of  this  program.  Map  scale  is  defined  by 
specifying  the  width  of  the  desired  map  in  inches.  The  maximum  allowed  width 
is  12.7  inches.  The  routine  computes  the  maximum  number  of  print  positions 
it  can  assign  to  each  grid  across  the  page,  and  the  number  of  lines  of  type 
per  grid  down  the  page.  This  is  accomplished  by  computing  the  following: 
COL  =  WIDTH  *  10.0  +  2.0 
where      COL  is  an  integer  number  and  equals  the  number  of 

characters  allowable  since  there  are  ten  characters/inch. 
PS  =  (COL  -  C0LS)/(C0L3-l) 
where  COLS  =  number  of  columns  in  matrix. 

FS  now  is  the  number  of  characters  between  columns  (see  Figure  h) , 
The  final  map  width  will  thus  be  COLMAX  =  PS*  (COLS-l)  +  COLS. 
It  will  never  be  greater  than  the  requested  size.  The  number 
of  characters  per  grid  is  thus 

XSTOP  =  PS  +  1 
the  number  of  rows  per  grid  (at  6  lines  of  type  per  inch) 
is  thuc 

YSTOP  -  ((PS+l)*6+5)/lO 
All  the  above  calculations  are  done  in  integer  arithmetic. 
Since  the  number  of  characters  per  inch  (ten)  is  not  equal  to  the  number 
of  lines  per  inch  (six),  distortions  of  the  map  result  except  at  scales 
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for  even  intervals 
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FIGURE      3.       FLOWCHART      FOR     SUBROUTINE     PRTC0N4. 
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FIGURE   4.      PARAMETERS     USED     IN     DEVELOPING 
A      PRINTER     MAP 
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where  the  grid  spacings  are  some  multiple  of  one  half  inch.   Table  1  shows 
the  distortions  that  will  occur  for  all  scales  from  a  0.2  inch  grid  (the 
smallest  allowed)  up  to  a  one  inch  grid. 

The  actual  map  is  computed  and  printed  line  by  line,  to  save  storage 
requirements.  A  line  of  blanks  and  plus  signs  is  first  generated  (STEM)  to 
define  the  column  spacing  along  the  top  of  the  map. 

The  routine  computes  the  estimated  elevation,  ZEST,  at  each  character 
position  by  linear  interpolation  within  the  lattice.  This  is  accomplished 
within  four  rested  DO-Loops,  as  follows 

(1)  Loop  through  map  matrix  row  by  row, 

(2)  Loop  through  current  grid  printer  line  by  printer  line, 

(3)  Loop  across  map  matrix  column  by  column,  within  this 
(k)     Loop  across  current  grid  character  by  character. 

Figure  k   shows  the  parameters  used  within  the  routine. 

Once  the  value  of  ZEST  has  been  obtained  the  routine  merely  has  to 
assign  an  appropriate  character  to  the  print  position  according  to  the 
contour  and  display  options  specified  (see  below)  and  then  move  on  to  the 
next  position. 

By  successively  estimating  elevations  over  an  entire  surface  for  a 
series  of  elemental  areas  (a  character),  the  routine  is  performing  numerical 
integration  of  the  surface.  For  example  by  subtracting  the  ZEST  values  from 
a  base  level  and  summing  the  differences,  a  measure  of  volume  between  the 
surface  and  a  base  elevation  could  be  computed.  This  is  not  presently  done. 
However,  the  routine  does  keep  track  of  the  number  of  characters  assigned 
to  each  contour  level  as  well  as  the  total  number  of  characters  within  the 
area  being  contoured.  A  frequency  table  showing  the  percent  area  within 
each  contour  interval  is  readily  produced.   These  data  are  further  analyzed 
by  Subroutine  HYPSO. 
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TAKLE  1 
MAP  DISTORTIONS   PRODUCED  ON  COMPUTER 
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Subroutine  FRTCONU  includes  many  options  allowing  for  various  types  of 
contour  intervals  and  various  types  of  displays, 
i)  Contour  Options: 

The  number  of  contour  intervals  can  be  specified.  The  maximum  number 
allowed  is  19;  the  default  option  is  5.  The  contour  intervals  can  be 
regularly  spaced  or  they  can  be  irregular.   Irregularly  spaced  contours  are 
useful  if  the  data  are  mostly  within  a  fairly  narrow  range,  but  there  is 
one  very  high  peak  or  deep  hole.  The  contours  can  be  specified  so  that 
they  are  close  together  in  the  common  data  range  and  widely  spaced  around 
the  peak  or  hole.  Regularly-spaced  contours  can  be  fitted  within  a  specified 
range  of  elevations,  or  the  range  can  be  computed  by  the  subroutine  to 
conform  to  the  total  range  of  the  data.  Statistical  contours,  with  «n 
interval  of  one-half  standard  deviation  and  ranging  from  minus-three  to 
plus-three  standard  deviations  may  also  be  requested, 
ii)  Display  Options: 

The  routine  can  display  any  type  of  contours  in  any  one  of  three  ways. 
A  series  of  c ont our  lines  may  be  printed  in  which  the  map  area  is  blank 
except  for  strings  of  characters  showing  the  location  of  contour  lines. 
A  series  of  contour  bands  may  be  requested.   In  this  case  the  map  will 
be  composed  of  alternating  blank  and  printed  intervals  since  all  even 
characters  (0,  2,  k,   etc)  are  printed  while  the  odd  characters  (l,  3>  5,  etc)  are 
left  blank.  Finally,  all  print  positions  may  be  requested  in  which  all  the 
characters  are  printed  without  any  blanks.  Figure  5  shows  an  example  of  the 
three  displays  produced  for  a  single  area, 
iii)  Contouring  Irregular  Areas: 

While  the  routine  requires  rectangular  elevation  matrices,  irregular 
areas  can  be  defined  by  setting  all  lattice  points  falling  outside  the  desired 


FIGURE    5.      MAPS   SHOWING   THE   THREE   PRTC0N4   DISPLAY   OPTIONS. 
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area  to  arbitrarily  very  high  or  very  low  values.  These  points  must  lie 
at  least  one  data  range  above  the  highest  point  or  one  data  range  below  the 
lowest  point.  They  will  then  be  classified  as  "no  data".  Figure  6  shows 
the  ranges  as  defined  within  Subroutine  PRTCONU.   The  frequency  table  for 
irregular  areas  shows  the  percent  area  of  the  irregular  area  falling  within 
each  contour  interval.  Figure  7  shows  the  contouring  of  a  drainage  basin 
by  PRTCONU. 
d)  Subroutine  HYPSO 

Figure  8  shows  the  hypsometric  curve  produced  on  the  printer  by  Sub- 
routine HYPSO.  Hypsometric  analysis  is  used  be  geomorphologists  to  measure 
the  current  stage  in  the  erosion  cycle  for  any  drainage  basin  being 
analyzed  (5,  8).  The  analysis  is  concerned  with  the  relationship  between 
altitude  and  areal  extent.  The  hypsometric  curve  is  merely  a  plot  of 
relative  altitude  (on  the  vertical  axis)  versus  relative  area  (on  the 
horizontal  axis).  The  relative  altitude  is  computed  as: 
Rel.  Alt.  =  (Alt.  -  A  )/(A  -  A  )  =  h/H 
where  Rel.  Alt.  =  relative  altitude 
Alt.  =  measured  altitude 

A  =  altitude  of  the  lowest  point  in  the  basin 
A  =  altitude  of  the  highest  point  in  the  basin 
Thus  the  value  of  relative  altitude  ranges  from  zero  (for  lowest  point) 
to  one  (for  highest  point). 

The  relative  area  is  computed  as : 
Rel.  Area  =  a/A. 
where  Rel.  Area  =  relative  area 

a  =  the  projected  horizontal  area  of  the  region 
containing  all  the  land  lying  above  the 
corresponding  altitude  within  the  basin 
A  =  the  projected  horizontal  area  of  the  entire  basin 
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FIGURE    6.     DEFINITION     OF    CONTOUR    RANGE    SPECIFICATIONS 
UTILIZED     BY     PRTC0N4  . 


FIGURE    7.      MAP   OF   DRAINAGE   BASIN    PRODUCED    BY    PRTC0N4. 
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FIGURE    8.      HYPSOMETRIC    CURVE    PRODUCED    BY    HYPSO. 
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Figure  9  shows  these  geometric  relationships.   It  should  be  noted  that  the 
area  below  the  curve,  called  the  hypsometric  integral,  is  a  measure  of  the 
volume  below  the  present  surface  relative  to  the  total  volume  of  the  solid 
formed  by  the  basin  outline  and  horizontal  planes  passing  through  the  highest 
and  lowest  points . 
e)  Program  CONTUR 

Program  CONTUR  is  a  modified  version  of  a  FORTRAN  II  program  originally 
developed  by  F.  J.  Rens  at  the  University  of  Michigan,  Ann  Arbor. 

The  main  program  (CONTUR)  utilizes  subroutines  SCAN,  TRACE,  CALC,  and 
DRAFT  as  well  as  the  CALCOMP  routine  PLOT  (with  entry  points  PLOTS,  SYMBOL, 
NUMBER,  FACTOR,  WHERE,  and  LINE).  A  dummy  subroutine  PROJ  is  included  so 
that  new  options  can  be  added  easily.  As  shown  in  Figure  1  Subroutines  SCAN, 
TRACE  and  CALC  form  one  sequence,  while  DRAFT  works  more  independently.   The 
first  three  routines  are  concerned  with  the  location  of  the  contour  line 
coordinates.  Subroutine  DRAFT  issues  the  actual  plot  commands  in  conjunction 
with  the  CALCOMP  routines. 

PROGRAM  OPERATING  CHARACTERISTICS 
a)  Program  MAPPER 

Program  MAPPER  will  compile,  load,  and  execute  about  50K  (octal)  central 
memory  locations  without  Subroutine  GRID  (i.e.  when  regular  data  is  being 
contoured).  When  GRID  is  needed  almost  70K  locations  are  required.  Running 
times  vary  depending  on  the  type  of  data  (gridded  or  irregular)  and  the  size 
of  the  grid  and  the  original  data  when  gridding  is  attempted.  With  regular 
(gridded)  data  an  allowance  of  2.5-3  seconds  per  map  plus  a  3  second 
compilation  time  will  be  adequate.  When  girdding  is  required,  times  will 
range  from  just  under  20  seconds  for  30  observations  being  interpolated 
onto  a  23  x  33  grid  to  over  3  minutes  for  1000  observations  onto  a  25  x  35 
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FIGURE    9.    GEOMETRIC     RELATIONSHIPS    USED    IN    HYPSOMETRIC 
ANALYSIS. 


C-25 

grid.  Figures  10  and  11  show  typical  deck  setups  for  gridded  data  from 
either  cards  or  tape.  Irregular  data  deck  setups  are  similar. 
b)  Program  CONTUR 

Program  CONTUR  will  compile  load,  and  execute  in  60K  (octal)  central 
memory  locations.  Execution  requires  about  i+OK  locations. 

Running  times  vary  with  the  size  of  the  data  matrix  and  the  number  of 
contour  lines  within  a  map.  However  for  commonly  encountered  maps  around 
30  seconds  is  generally  adequate.  Figure  12  shows  a  typical  output  from 
Program  CONTUR. 

PROGRAM  LISTINGS  AND  DOCUMENTATION 

The  following  pages  contain  listings  of  the  various  routines  in  order. 
All  are  written  in  basic  FORTRAN  IV.  Each  routine  is  documented  through  the 
use  of  appropriate  comment  cards  at  the  beginning  and  throughout  the  program. 

Each  program  has  been  analyzed  by  Program  TIDY,  available  as  a  library 
program  at  Purdue,  which  resequences  all  statement  numbers,  removes  extraneous 
numbers,  and  generally  makes  the  program  neater.  For  this  reason  the  statement 
numbers  in  each  routine  increase  progressively  in  increments  of  five,  and  all 
FORMATS  (if  any)  are  collected  at  the  end.  The  cards  are  each  given  a 
sequence  number  which  appears  at  the  far  right. 
a)  The  Use  of  FORSTAT  Data 

Each  routine  was  next  analyzed  by  another  program,  FORSTAT.  This 
checks  the  syntax  of  source  (FORTRAN)  programs  and  produces  cross-reference 
listings  of  all  variables  and  statement  numbers  found  within  the  routine. 
The  cross-reference  listings  utilize  a  series  of  "FORSTAT  Numbers"  which  are 
found  at  the  far  left  of  the  listing.  A  vertical  ink  margin  line  has  been 
added  to  separate  these  numbers  from  the  rest  of  the  listing.  The  cross- 
reference  listings  are  found  at  the  end  of  each  routine. 
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FIGURE    10.  TYPICAL    CARD   DECK    SETUP   FOR    ONE    ANALYSIS 
USING    CARD    DATA 
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FIGURE  II.   TYPICAL  CARD  DECK  SETUP  FOR  ONE  ANALYSIS 
USING  DATA  STORED  ON  TAPE3. 
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FIGURE  12.    SAMPLE  OUTPUT-  PROGRAM  CONTUR. 
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All  symbols  defined  in  the  routine  are  listed  alphabetically  gnd  are 
identified  according  to  type,  definition  and  references  (where  referenced). 
The  following  symbols  are  utilized  in  the  cross-reference  listings: 
Types  I  =  INTEGER 

E  =  REAL 

D  =  DOUBLE  PRECISION 
C  =  COMPLEX 
L  =  LOGICAL 
N  =  NAME  LIST 
ASF  =  ARITHMETIC  STATEMENT  FUNCTION 
Definition  C  =  COMMON 

D  ■  DIMEIBION 
E  =  EQUIVALENCE 
All  other  terms  are  believed  self-explanatory.      These  cross-reference 
listings  will  be  most  useful  should  any  reader  desire  to  make  changes  to 
these  routines. 
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b)  LISTING    OF    PROGRAM    MAPPER 
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ImIS    PWOi.fA"    cOnIaInS    In«tE    SeUIONS.oNf     IS    USfcD    *HrN    IhhfjGULah 
iia'a    IS    SUPPL  ItOtiNOIHtH    ""En    OHIUUtii    Ua  I  a    IS    SUPPLlEUi     aNo    IMF 
TmIHIj    CALLS    SUbH'H'TlNt    PhICOn»    TO    P«Or)UCF     IMF     MAPS. 

1p*FGULaw    [)a1a    m«Y    Ht    SUPPLItU    FROM    F I ThEW    CaPUS    CM    IaPciTapF?) 
IT    lb    GAtUOfcn    hV    SLiHnuuI  [ML    ii«lD    aCCOmuImi    In    SPFcIf  ICaT  IOnS    SUPPII- 
Fl)    oY     TMt    USFP.     Slf'Ct     SUriWUUTlN£     bHlu    CANNOT     HANOLE    MOPE     lhAN     1000 
ohSLhVaTIUns    aI    a    TlKti    phUOnaM    maPPE"    GE ogpaPHIcaLLY    *P*E-SopTS* 
I  APbtH    APHAYS    IMU    SUbSE  I     IHtPPEStNl  |w.    SlMAPLAS)     tt»Ch    COnTaInInG 
I.FSS    Than    1000    PnII»TS.     TmE    COMPUTED    GHIOnED    VALUFS    MAY    HE    PLACED 

ON  A  IAPK   (lAPt3>   ("H  PuNChtU. 

GHluutO  UAlA  E'AY  »E  SUPPLIED  FHO"  EITHER  CAPOS  OH  TaPF(TaPE3). 
SINLt  THIS  DATA  Akt  (.LhtAuY  IN  A  SUllAhLE  FOHM  FOR  SUBROUTINE. 
PWTCONI,  ,  NO  FUKTHeH  MOllIE  ICa  '  10^5  A«E  Hf'JUIPtnt 


nfSCHIPTIUN  OF  conTKOL  CahuS  — 

I)  U*,TA  IDENTIFICATION  Cahu  - 

iHf  FIRST  CONTMOL  Ca«U  IDENTIFIES  Thf  TYPE  OF  OAlA  mFINO 
SUPPLIEn  nY  FlThtH  OF  l«u  CODE  ».OPL)S  IN  COLS  1-10  . 
•IHPeGUlaP  -  lwDlCAltb  IHH£r,ULAPLY  SPaCFO  OaTa  . 
•G-tluDtD    -  UDiCAItS  iiHloiitD  DATA  . 


A)      LAHOS     EOH     lKhK.Uj  »►     l>AlA 


Al)     MaSIEH    OPIIOn    LAHLi    - 
nTapE       (COLS       1-    b» 


KO-ilUP     (CCiLS 

ioji 


nOhS 

NPT 


Int  »OHO  *IaPE2»  IN  IhESF  COLUMNS  INDICATES 

InAT  IMF  DATA  IS  TO  MF  SUPPLltO  FPOm  TApF. 
1U )  -  fnt  »0Ho  »Samf*  IN  TntSF  COLOPN5  iNnlCAlFS 

IhaI  IMF  KtllUESTEo  "*PS  a«E  'C  ~,t     IutNTJCAL 
(COLS  H-lb)-lnt  »0Hn  <UPEJ*  In  ImEsf  COLUMNS  iMiICATtS 

InAl  thf  COmpiiTei)  GPIU  VALUES  APE  To  h£ 

PLACtD  ON  (APE   , 
(COLS  )6-i0»  -NUMHtH  OE  SUHAHFAS  DATA  IS  TO  MF  PHf.SOPTEO 

iNlo.  TmIS  MUST  *E  aI  LFaSI  ONE  EVtN  IF  NO 

SUBDIVISION  IS  PFOOIHFU. 
(COLS  ?l-2b)-N0M,tH  of  O.iSFHvaT  InNS  hfING  SUPPI  Ihi 
(COLS  26-30)  -nUM.tH  of  OhSF.hvaTIonS  HFCOHOE  0/CaHii. 


10 
12 

I  H 
lb 

1M 
20 

?i 
?- 

?o 
21 
30 
32 
34 
36 
3« 


So 

S2 

s. 


Hi 


106 

l  Oh 

10 

12 
1« 

In 
llH 

l?0 

l?Z 

I2A 

1  ^n 
l?d 

130 
132 


TMF  fOLLO«|NG  CAWDS  HE.OUIMtO  ONLY  IF  OaTa  SUPPLIFO  EOm  C«hOS 


C— 31 


l/)     FOPmaT     CAUD    - 
FOrtM  (COLS 


1-/0)  -tUH'Ml     FOP    DATA     CAPIiS, 


a  j)    akea   hf.adfp   r/iHU  - 

fac'uh    (cols     l-ju)  -AL^MUNUMtmc    iufnT  if  icat  ion  ok    facToh, 

A«Ea  (COLS     A  WO)-ALP*ANU"tHlC      I  Of  ~ll  1  F  1  C  A  T  I  ON     OF      APFA. 

I  f  r  (--t       (      OH  s    to    1-1F    maIa    IS    IoPOoHaPhIc    ofS  UjuaLS  .  I  I  ype«I     . 

ik   i >- 1 a  MtAsuPts  rosis   .    i  r TF-t  ■<;   . 

IF     ijaIa    *t«SUPES    HE"tFITS» 1 IYFE«3    . 

CA)      OAlA     CAPOS     « 

THE    FOLL't»lN(i    (  A"OS    A*t    UbtU     10    DEFINE.     T  -it    SumapeaS    - 

aT    LtAbl    UF.F     St  I     IS    aL«AYS    HtiiUI^tO    i    F«F'I    »Mt»l    NO    Sum,  I  vision    is 

f.ECtSSA^Y      .      TnF      OVFH«LL     uMtA     MUST     Al     I. FAST     HF      DEFlNtU     CN     THtSF     CAPOS 

AS)      SUHAmEA     IIILF     LAoU    - 

IIU  (COLS        1-/UI-aLPmaNUmEP1C     IlJF  NTIFYINfl     TIILE. 

Ah)      SUHA«tA  LI"II  SPttlF ICAl IUN     CAHU     - 

1M««  (COLS  1-1U) -MAA|KUM    A-L00PU1NaTf  OF  S0n»HfA. 

Xm]n  (COLS  1  1-iiC) -MINIMUM    A-COOKUI'JaTF  OF  SUrtA*tA. 

Yma*  (COLS  ?1 -JO) -HA* I  MUM    Y-COOPOINA'F  OF  SUBA«EA. 

YMJN  (COLS  31-*0)-MlNlMUM    Y-COOPDINaIf  OF  SU1AHtA. 

87)     UPlU    SPECIFICATION    CaHU    - 


tlUO* 
NCOI 

(.Slit 
MM) 


|PUN 

NMAH 


(COLS         1-J     (-NUMHtP     OF     «U*S      IN     GPIO. 

(COLS      »-6    )-numhEF<    OF    COl  UMNS    If,    «m Tfi . 

(COLS       7-lO)-siZf;    OF    bWlo    ]N    »-Y     InPut    OA I A    UNITS. 

(    C  01        11       )-lF    nNO»0»&K|f)    mOUnO/ihIES    CALCLLaTfu    FHOm 
UAlo     ACCUSING     TO    Al  OOFtJ  TF>M, 
IF     uND«lt(iF<ln    M'iONOAr(IES»«rtA*f  XMlNt  YMAAt 
TMll,    AS     SPt-CIFIFO    PPtVlnilSLY. 

(       C"t     U       l-IF     ILIST»1«)NPUI     (IPheouLAH)     CAlA    LISTED. 
IF     ILIST"»2«ouTPuT  (GPlOuEol     DATA    LISTED. 
IF     ILIST«3»H0TM    InPdI    AnO    OUlPUI    DATA    APF 
LlSltu. 

(         COL      13         )-lF      IPUN.l.UWIO     VALUES     APE     PUNCMFO. 

(COLS  1»-1S)  -nUmmlm  nE  map  COPIES  TO  HF  PPCUuCEl)  FOP 
l»lS  SUMAHtA, 


Cumplf  IE  COnTmol  CaPU  PACKAGE  «Ph  CAwdS  in  Pa«I  C. 


u)  LaPDS  FOM  GPIulitO  UATa  -- 

u)  I  maSIF'H  OPT  JO'   CaHU  - 

LIST    (   COL  1    >-U  LlST.l.  jnpuI  oaTa  AOUAY  IS  LISTfd. 
"(map    (COLS   2*3  l-nuiiHtK  OF  map  Copies  To  Pt  PPCOUCEU  F  OK 

InIS  DATA  St  T. 
«')NlUh  (COLS   6-l0)-'"t  *UMM  *Sa"E*  In  ImeSf  colupns  indicates 

ImAI  IMF  PtOUESTEU  "APS  APt  10  ,,£  lutNTJCAl 
maTw|«  (COLS  ii-1S)-ImE  «0M|>  *1aPE3«  IN  IHESE  COLUMNS  INDICATES 

TnA'  IMF  UAlA  IS  TO  nE  SUPPLIED  FP0«  TAPE. 

ThF  FnLLO*lNr.  CAMjS  KtuUIPtU  only  je  oaTa  SUPPLIED  F0M  CaPdS  - 


'7)    fdwmaI  Capo  - 
E08m    (COLS 


1-/0) -F UW~»T     FOP    DATA    CAPI'S. 


i)|     UATA    SET     hEAiiFF.    CAhU    - 

UU  (COLS       l-7U)-«LPMAivUMtPlC    lOFNlIfYlNfl    TITLF. 

'jC'JL  (COLS     71-73) -IvUMHtP    qE     COLUMNS     IF'     THE     AHHAY. 

nP0«  (COLS    7*- lb) -NUMHtH    OF    KomS    IN    Tnt    ABPAY. 

ITYPF        (        CIS    7H     )-|F     0»IA     IS     I  OPOr,PAPH  1  C    uf  S  I  DUALS  .  I  I  YPE  «  1 

IF    ij«H    FEASUPEs    COSlS    •     |lYPE*i?    . 

IF    ii.IA    "EASUPts    HE' tEITS.IIYPE.3    , 

LUMPLFlE    CCNTpOI     CAPO    Package    *IIm    CAPOS    7'«    PA"T    C« 


1  M 
I7u 

i  r? 

i  r* 
i  ih 

17H 
IHU 
IflZ 
18* 

1  ib 

1  HI 

19U 
kV2 
19« 
I  wn 
lv- 
20  U 
202 
20<t 
20A 
20h 
210 
/\l 
21« 
216 

ZlFJ 

,vo 

Z?k 
226 

<>?* 
21g 
B\? 
23« 
^?A 
?1H 
<!mU 
2»i 
244 
<!*h 
2*« 

<"SU 
?Sl 

?s» 
?=>>, 

?SM 

260 
264 

^F.6 

2iSH 
2T0 
211 

?7«. 
?Th 
?1H 

2H0 

?H^ 
?>)<. 

?«« 

?HH 
290 

^92 
?9« 
296 

^v« 
300 

?0? 
JO* 
30S 
On 

no 
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spf  f  if  l  C/.  i  I uk   tA-io    — 


n*aH    ^Ms    OF     Thesi     wt'-ul-tu    UNLESS    *SamE<    FNTtHEn    ON     If-E    aPPHuPPUU 

rmuniA  of    c"MW"i    c^us   <-i    ri«   hi    . 

T|)      »«K     sftCIMCM  ION     LAWK     - 

f  II  J  (  COL  1  -^  )  -  T  Mfc      uOMntx     OF      CONlOUH      IrTfHVALS      |MA«lPU"      =      Uli 

I    (-yn    =    l.aIu    HAi'Ot    /    UEMHEti    COnTOUH    InIFRVal     I. 

lU'jHiCOLII-      IF       '"OH     =     U>     F>HMfiHAM     USFS     Mi*|MIIM     ANU     MINIMUM     FM)M 
DATA     FU.<     UPPtH     AMU     LOWEW     C<1K  I  OUH     L  I  M  M  S  . 
IF       lUUK     =      1,      UESlMtll      TUP      AMI-     H'lMOM     LCNlOUHS     HEAD 

It.    EHOM    SefcXiAL    PHUiWAM    CA«r      (SFF     BELOW). 

IF      IUUH    =    <;,     VahIABIF     CONTOIIHS    "F»U     in     ilo«     10 

HIGH).        I'*     I"1S    CASr     SEI     COf     EUUAL     (0     T-E     NomhEh 

of      »«LUti     r  ■  I     HF     HEAD     IN. 

IF      IHOH     =      I,     COnlDUH     INIFHVaL    tlFCOMtS     1/2     SIanoapu 

I  >F  V  1  u  f  1  I  !.>      I      "AX     =     ♦(     SIl,«Al      "IN     =     -J     5I(,«A      ). 
LINfSICOL*)      II      Ll'*ts     =     »•     CONlOUH     LINfS     AKt     PMnIfO. 

IF      I  |l«tS     -      1  .      Al    lEHI'ATfc       (EVEN     VAlUtOI      C'lNJUoH     BANDS 

AHfc      HMlNltU 

IF     LINES     =     ■><     AIL     OA'Hib    PMjFlED     (ALL     PRINT     PnSNS). 
wiiMicnL    "!-•*)-    *IuT"    is    mam    wliiih    In    UGhEs.     If    GHtATth    Than 

1  2  .  I    UH    nLA'.M    st  I     1(1    10,0 
InYPSulOU      I  <  >  >  —  It-     lnVHSu    =     I  iSUmHUUI  I  Nf     HYPSn    CALLEn. 

V?\      SPECIAL     f«nis    -     HtuUIKtU    UNL"     »"l»l     ToUH«     1     OH    2     • 
IF      IOUH.     =      I 

flNF      CA"H     C^Nl  A  lNJNii- 

/<-  Ih  (LOLS   1-101-mj.nimUm  CONTOUH  flc.vaTj.un. 
/.<AX   (LOLS  11-iO) -MAXIMUM  CoNTuUH  ELEVATION, 
IF   IUUH  =  2 

ONF  C«HO  Cu"l  AtMNli- 

INCON     (LU!  S     1-F>0  > -F  OmmaT    kik    KULLOWINC,    CaTA    CahdS. 
PLUS    0'  F      UK    Hunt     C<."US    LONTAINING- 

A  stoltb  OF  VuLUES  HEHkESf  NT  Uo)  Thf  MElC-HlS  OF  THF 
l«i'|b"LAHi  T  ShmlEo  COnTOUhS.  TnsE  MUSI  EE  apHanuEii 
I  uuh    LU»ESl      III    hIr,HF.ST     IN    OHUFH, 


lOTt    -- 

TmIS    PHObhA-    UTILIZES    LaoELLLEU    COMMON 
/  r,  1 1  a  T  A  / 

/r-APS/ 


I'EUUlHFIJ     SUHH0UT1OFS    " 

(•H[|)     -     (1NL  Y     *F<E"'      IHHtUULuH     |,AlA     SUPPLIEO. 

PH  TC0N4 
STAIS 

OPTIONAL   SUHHOUTINES  - 

HYPSO 


OlMFNSlUN     AHFAIil,     UCTU«H) 

DIMENSION     TIILI7).     FUHflM.     ZlSO.bol,      ICOIJFlS),     XTllUOOll     T  T  (  1  OtlO  ) 

li    /tuouOIi    «l|(KiO)»    yiiuuiii.    wiiouni 

INI  f  f,F  H     CON,  |  OljH 

[HlKGF.H  SYSlN.srsol'SOF  I  tS()E2»SUF  J 

common  /GDaTa/  koon  I  .  *,  r ,«,  I  [Tl  >nf<Ow,ncol»  I  I  YPf  ,6S  I  ll  ihNU<  TL  I  ST  ,  I  C 

IUN.K   lAPt,XMA»,XMIN,YMA*,TM|N 
CO  "MOM     /maps/     rohllIOUH.LlNFS,w|0TM,IHYPS0»TNC0N(F')i0PAN(i(<!u'),ZMlN« 

1/MAX 

oaIa    SYSIN.SYSnT.SUE  1  •  SuE <!«bUF3/5ib<  1  *?tj/ 

DAlA      (  ILOHF   I  1  )  «  I  ■=!  ib) /SMSA"t     •  S  Fi  T  APF  3  I  I  I|M»  I  HHE  Tl  IL  A  H  I  1  0"  *C'H  1  OOE I '        i 
ISHl  APF  /•/ 

FuulVALt^Cf     (/mTli     1/(1001 1  iYT>  i     '  /  I  200  1  '  <  1 1 ' 

HE"  I  NO    SI  F  J 

HE».|N'l     SUF? 

N  I  AHF  all 

K  Ul    F  =  ( 

HEAli     (S»S1N.??H)      IS»LM 


M« 
31" 
11« 
3?0 
l?d 
3?* 
I?h 
I2rt 
33'l 
33(- 
334 
33b 
33>« 
34  I 
3«^ 
IM 
3*b 

14- 

35U 
354 

ihh 

35H 

Ifcil 
3fe^ 

)h4 

if)H 

17" 
37^ 

3  ;- 
I  7ft 
I  rn 
)H0 
3H? 
3H4 
3FCJ 
3HH 

3s"u 
34i 
I  j» 
396 
34H 
»0') 
»02 
404 
4llh 
»0« 

410 
412 
A  1  4 

4  IF, 
41M 
420 
4i>2 
424 
»^^ 
4?>1 
4  30 
432 
4  )4 
4  11 
4  (x 
440 
44  2 
444 
o  w 
444 
4S0 
4S<> 
•  ^ 
4S6 
4SH 
4A  i 
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r  M  AU     CA-<I>      I  [)F  ■.  I  !  t  *  I  NO      I  Yi't     Ut-     DATA 

c         

If-       (F    >F  .SYS  IN)      1 75i 10 

in  Ir     (  ISwCH.Fu.Inil'E  ( J)  )     M)Ofc'  =  l 

IF     (  lS.vO.t-0,  K  ii|,fc  <»)  ]     rui ^  =^ 
(•0     To     (  15»  1  1  f.«  1  >•(■  )  •     imji   (■ 

(         

C  iHhU.'ii  ftw    DATA    Sk-(  I,  Ir  ILU 

c         . 

C  WFAL)     "  A  S  t  t_  h.     nPllilN     i>Pt  C  IF  It-A  I  1UN     CAKij 

r        

IS  Ht»M      (  S  r  S  |  r   ,/  )i    )      Ml  ,.PF.  ,i\iji»I  UH,  !  (IU  1  ,NAW|-  A  ,N()BS  IMP  1 

If       (I  Oil  I  .(■!,.  |  Cil't    <<->  )      X  I  API-  =  1 

IF       (MaPF  .'  >   .  Il  "I   F  IS)  I      >,<J      I  J     40 
c  

(  T.iPO  I     Hi  I  /i     F  ►.iir.      I  r.FF  r 

C  i-t*u    IaMFL     AMI    llll»     l<hCUwub    h"OM     I  AHF  ? 

C  AMI     H«rt      L'AlU      HllllHi.S     0  j      lAKt.1 

<  ---- ..... .. 

MLull      (Sjt?>      F  AlTdP.AKtA  t  1TYPF 

IF     uio'-s.i.r .  I  mil  i    (ju    fu  io 

Kivul=l 
r  lf=N,)„S 

h  I  s  *  =  1 

i>u    In    e-s 
r 
^i/         KNu~i=wjds/i  ono 

JUU":  gU'iS-K  Ml  'M«  1000 

K  I  -,  =  1  DUO 

nl'j.rll 
/S  llO     IS     ^  =  1,1.  MUM 

<0  "MM     ISHf^l      1*1  U  )  lY  I  (1)  >^T(1)  »1  =  1  <lvlFt) 

is  ..hitf    (SOFu     ( <  i  (  1 )  .r  r  (  1 1  ,/ i  (  i )  ,  ui  ,nt..  j 

le     HI  s».f.j.  j  )    (.o    To    7n 
r^  I  h  =  jnmjm 

K    IS»=1 

Cl>     111      I  (J 


t*  7  U 

hi* 
«  7h 
»  71 


SOU 

502 
So* 

506 
SUM 
SM 
512 
M» 
Sift 
Mi 
5?0 
522 

S>H 

5?6 
5?H 
Sll) 
S  )2 
5  <.» 
536 


ni-  1  A     F  sop    CAPOS 

"UU    Kim*a1     CfSO     Hifcl*    aHl  „     rtt«|)£P    CA"Ll 

Thfn    »un    i.«TB    F»"»     CAP-uS    AN.)    PL«Ct    ON    TflPf] 


5*0 

5** 
546 


Ht»ii    is'ST*  .^nn)    fUK« 

Htnl     (sYsfMlHS)     FAL  1oi.»«hf  ».  I  TYpt 
I  F      (NOMb.M  .  I  llm,  )     (.u     I  li     15 
JNUmsv,,)mS 
(,0     10    i" 

«.NU"  =  NO'IS/  |  00(1 

jNU'"  =  MlinS-|.NU"»inOU 

NblFP=(  1  000  ,/y  |  i,aT  IimpI  )  )  ♦.  v 

UU      SS      K  =  l ,<NUM 

!Sm,=  1 

HO    SO    N-  1  ,»M  I  P 

If  Wis  |hFi--*NP1-j 

wtAO     iSVMI   llUHH)      (Al  (  I  |  ,YT  (  1  I  ,Z1   1  I  I  ,  1»IhfG>  ItSOl 
|HF,i=  |  hf  r, »npT 

wKITt      I  SUM)       (  «T  (  1  )  .  VT  {  1)  ./  I  (  I|  »  1=1  ,  1000) 

N5IF.P=(FL0A1  IJi  iif'l/FLUul  <NPl  >  )  ».9 

lntr,  =  l 

nu    ss    '4  =  1  '"  S1F  p 

IF  V)=lHF  h»NP 1  -  1 

wFAn     (SYSlN.M  KM)      (<l  (  1)  .YT  (  I)  IZT  (  I)  .  I«IHEG,  IENOI 
lHt(i=llitbONP  I 

»hiif     (bliFl)     (XT  (I)  »YT(i>  t/l  (I)  •I  =  1.jnhm) 


siiMlMVlOF     I  "l  1  A     IiTh    5UHAl-t»S    oY     «;tAWCMlN(,    OATA     Ow     TaP£1 


H«U     MM«"M      rillF      CA«lnfHtn     SUH»KtA     LIMITS     TAPD. 

THtrvi    (,Mln     5PFC1F  1(  ATlON    {amij. 

FAIM     SIMnwtA     "I'Sl      I'AVF      Nil     MIlMt      THAN     1T00     OHSFhi/ATTONS 


SS" 

SS? 

SS4 
SSh 

SSM 

sso 
562 

SS4 

Sss 

SftM 

570 

S7^ 
S74 
S7h 
S78 

5«0 

SH? 

Srti. 
S"h 

590 

S<J? 

sv* 

596 

S9M 

sno 
602 

*  n» 
606 
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65. 
66. 
AT. 


"1  . 
st?. 
-  3. 
a*  . 
"5, 
u6. 


95. 


100, 

in  I  . 
lo2, 

In.). 


In  4. 
105. 
I  Oh, 
10'. 
I  OH, 


1>U      1  7(1     NN  =  1   ^4"f  A 


c 

id  a 
c 

c 
c 


KOU> 
Kf  <I 
«F  An 

Ht*i> 
HE  All 
IF      ( 


N(j 
I  5 
IS 
(S 

NOH 
=  1 

NOR 


<;nF  i 

YSIN.200)     T ITL 

YSIN.P35)      AMAA,»X»lN,  YMAA.  YMlN 

YSlN.P'-OI     nhuw,nCOL,GSI^F-.HN(>,ILTST,Ipun,nf<AP 

«;.(,!.  100U1     6'i     lo    7  b 


*TS«  = 

f,0     TO    HI 


KNOM« 
JNJ«  = 
K  1R=1 

KT5w« 

no  *b 
pt 

on 


NOHS/ 
NOHS- 
0  0(1 

0 

K=] 
AD     (5 

9n    i 

If    ( 

IF     I 

X  (K0 

Y  l«fl 
w  IRC 
KOUN 
IF      I 


1  0  0  0 
MiUF'»1000 


DUI      (AT  ( 1  )  iYT  (1)  ,i!T  (I)  .1  =  1  iKTK) 

=  I  ■  f .  I  b 

XT  ( ] )  .L  I  .XMlii.ow.xI  ( 1 )  .gT.XmAX)    GO    TO    90 

Y  11  I  )  .1  I  .YM1N.UR. YT I  1)  .OT . YMAX>     GO     TO    90 
UM  )  ix  I  I  1  ) 

lint  i  =y  r  1 1 ) 

uMi«;Mii 

TaROUM  «1 

KOUNT.GT.  1000)     00    10    105 


IF  PMFVIUUS  sTaTfxINT  IS  IhueiMORE  Than  1000  POINTS  HAVE 
hEEn  FOUnU  In  !hf  (  UhHtNT  SOhahEa.  ALL  (JKIOOT^G  ANO 
CONTOURING  MTtMPlb  A*E  AnOHTLU  FOR  THIS  AREA. 

hut  map  control  caros  ahf  KtAu  so  that  correct  control 
CaRUS  will  he  foiinu  hY  siihSEuuent  analyses. 


CONTINUE 
CON  I  1NUF 

IF  (RlSw.tO.l I  GO  To  100 
kTr»jnum 
kTSw=1 

r,n  ro  hs 
rOUn  (*kounT-  1 
CALL  SUHHOUTINf  r,Hlo  10  GHlu  IhE  OAIa 


CALL  GRID  I/. SO, SO) 
GO  10  IJ5 

c 

ins      myswsi 

WRITE  (SYS0T,^<,b)  T  J  I  u 

SO  TO  IIS 
C 
c    

C  iiRlUUtn     INPUT     [iATA     SptCIFltU 

C  

C  RFAO    MASltR    OPTION    SPECIFICATION    CARU 

C  

110  RFAO     (SYS1N.I90)     LIST >NmAR,kOnTuR,MaTRI X 

IF     ItOF.SYSIM)     175,11b 
lib  IF     (NPAM.LT.l)     nmaPiI 

IF      (MATRIX.  fcO.  ICUUE  (<;)  )     00     TO     1?5 
C 

L 

c 
c 


CaRU  inpuI  data  SPfrlFItu 

READ    FORMAl     CARD.     THtN     T I  1 Lt    CARD 
then    The    matrix    row    hY    row 

HEaii     ISYSlNi^On)     Fowm 

RFAI)     ISYSIN,  195)      T  I  TL.NROw.NCOL.  I  TYPE 

1)0     120     I«l,r>ROW 

RFAI)     (SYSIN.FORM)      (  2  (  1  ,  J  )  ,  J.  1  ,  NCOL  I 

GO    TO     US 


1NR01     DATA    ON     1ARF3 


uf.AU    IiiEnT  if  icatiof,   hecohu 

TMFN     RFAO     MATHIX     ROW    HY     RO. 


MO 
612 
M» 
616 
lib 
620 
622 
624 
6?6 
628 
b)« 
6  32 

6  14 
63b 
63B 
b4  ] 
6*2 
644 
646 
648 
650 
652 
654 
656 
658 
660 

bb^ 

666 
*,-,§ 
670 
672 
674 
676 

67h 

680 

bHC 
6h- 
6H6 
688 

6<J0 

694 

6SI6 
f.g« 

7  0  ij 
702 
70* 
7  0  6 
708 
Mn 
71? 
714 

716 

718 
720 

7?? 
7?4 
7?6 
728 
730 
732 

7  34 
716 

738 

740 

74? 

744 
746 
74H 
750 
7S? 
7S4 
7S6 
7SH 
7  6  0 
76? 
764 
766 
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109. 
110, 

111. 


1)2. 

113. 

1  1*. 

11*5. 

Lib. 

11?. 
MB. 

119. 
12(1. 
121. 


122. 

123. 
12*. 
I2S. 


126. 
1?'. 
12*. 

129. 

|3U« 

111. 
132. 

L33. 
134. 
13b. 
136, 
137. 
13H, 
139. 


1*0, 

141  . 

1*2. 
1*3. 
1**. 
l*b. 

14ft. 


1  JO 

c 
c 


i'MIi     (SnFJ)      IT  I  IL  1  I  I  •  1  =  1  <4)  ,  I  Tl  II  (  1  1  ,  I*>m7)  ,NCUL  ,l>BO«,  |  TYPE 

IJ(i     |  ifi     |  =  1  .1  HOW 

Kf  Ah     (SlIE  jl      (/  (  1  .  J)  ,  J=1,"(.0|   ) 


prfPukE    cunTOIik    "HIS 

UF  AO    CONlOUB    map    CONIHULS 

TF     *S»"t»    OPTION    xf  UUhSTt-Ln    o»LY    UNt    SET    OF    tONPOL    CahdS    IS    NEAO 


1 3b  00    \bb    N  =  ] «f "AP 

If      (KOl-.tUW.E*.  ICnUL  (  1  I  .ANU.N.OT.)  )     sr     To     1*0 
»tAI)     (SYSIN, 21)51     CON,  IOUW, LINES, «|(1TH,  IMYPSO 
IF    ( TnuP.M  ,«!.un,  iuum,le.o>    oo    To    isn 
r,o    Tn    Il*o,  l*b)  ,    Iuum 

141)  MEAD     (SYSIN.210)     /Ml-WMAA 

(,u    to    lin 
c 

14b  DtAO     ISYSJN,  <JUU)      INCorv 

I'tAO     (SYSIN«|NCCH)      (O-ANOH)  .1  =  1, CON) 
ISO  CUM  INOE 

c       — - — 

C          It-     LIST     i     1,1   1ST     H.PUl     LOilA    MaTHjik 
c  

IF      (N.f.T.l  .OH.LISl.Nf.il     GO     li)     1*11 
«"1 Tf      I Srsul ,21b)     I  1  It 

llll      ]SS      Ir  1  ,NMU» 

1SS  «-UTf     (SrSi'T.220)     U  ( 1 ,  J)  •  J*l  ,NCiiL) 

c  

C  ruLL    PP1C0\4     To    pPO[}IJLt    (.UNlnOH    MAPS 

C  nmap    SPECIFIES    N"«ht"    OF    '"«F>s    Id    nF    PI'OUOCED    f  OR    THIS    LAlA 

(  . ... 

Ihfl  CALL    PHTrOTi*     l/,T  1  IL,'ihUwiNCOLi6OiS0) 

IbS  Con  I  1MJF 

I  In  C^T  INut 

HO    TO    b 
C 
i  fs 

1  Mil 


1  us 

I  VII 

Ivb 
200 

F05 
2111 

/  I  5 


220 
22b 

23u 

Z  <S 

/4    i 

/4b 


Ht»  INO     SOF  i 

COM  I  Not 

S  lUP 


FOhn 

FOHk 

FUH* 
F  Ui- 
FOM" 
F0«> 
lHt     ( 

2,29r 

3F  OH 

4HASF 
F  UK" 


F  UH» 
F  U«" 
FOH> 
110,., 

END 


T      [  jA)  II,  1  0«,3AlO,hA,  I?  I 
T      III  .T2«2*»2At>) 

i    ( /  a  i  n ,  a  1 3  •  i  i  i 
I    I  Ml  II) 

I      (I?,?I|.F<..u«I2> 

r   (2Fio.O) 

T     I)H1 ,5S*i29h ./,StX,29HI  T 

AwS    SVSIF^  I  »/»3.<>X>29Hl     PUHOOE     UNIVeHSJTT     VERSION     I./.S6* 

--i///, *S«,30HUFvELOP    CONTOuH    MAPS 

■MuuEn    UATA.  t///,?b*»  ibhTHlS    map    Sho*S    »7»l0»"i*0*«3bHlT    IS 

Of»    1HF     FOLLOmlNO    LIxlA    --,///) 

r    i  iho,  i?f  i ii.i,/ 1 ia, i/f  10,3) ) 

T     (Alii) 

I      IJAS.31SI 

T     UF  10.0) 

T     (21  1,F4.n,3l 1 ,12) 

T     I IhI ,4^M•»•FHH0H*••     100    many    POINTS    IN    TEST    AhEA       >/,12(«7A 

12A,2SH(,OTNl.    ON     TO    fiF_«l     PfOOtST.) 


7  7u 

ru 

77» 

lit, 

77m 
71') 
7H2 
7m« 
7^ 
Mn 
790 
792 
79* 
79ft 
79d 
«nu 

F302 
RQA 

H0»J 
HOH 
H  I  'J 
Ml2 
HI* 
H16 
Ft]H 
H20 

h?2 

H2* 
H?« 

H2n 

M   ill 

H  )Z 

H  1.. 
h  16 
H  1H 
Ft*0 

M4^ 
R44 

M4fj 
Mi,  M 
hSii 

eiS2 

MS* 
MSft 

Mb« 
660 

Mh^ 
MM. 

Mf>n 
m  7u 

H7<> 

M   /«. 

H7f> 
H7n 
MH'J 
6H2- 
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c)     FORSTAT     CROSS-REFERENCE    LISTING    FOR     PROGRAM    MAPPER 


PROGRAM      MARRfR 
STATEMENT  NUMaER    OEFINIMON 


REFERENCES 


5 

1  3. 

129. 

in 

17. 

lb. 

is 

20, 

19. 

2o 

29. 

2*. 

?s 

33. 

28. 

30 

3*. 

39. 

3S 

35. 

33. 

♦  0 

40. 

22. 

*s 

45. 

42. 

so 

53. 

50. 

ss 

5*. 

4b. 

60 

ss. 

44. 

bS 

60. 

57. 

70 

62. 

3b. 

7S 

74. 

69. 

80 

78. 

73. 

MS 

79. 

93. 

Sin 

HH. 

80. 

81. 

9S 

h<v. 

7b. 

ion 

9*. 

90. 

los 

97. 

87. 

lln 

100. 

19. 

lis 

102. 

101. 

120 

107. 

1U6. 

125 

109. 

103. 

130 

111. 

110. 

13S 

11?. 

9b. 

99. 

140 

117. 

116. 

1*5 

1  IV. 

Mb. 

ISO 

121. 

lib. 

118. 

lbs 

125. 

124. 

160 

126. 

113. 

122. 

165 

127. 

1  12. 

1  70 

128. 

62, 

1  7S 

130. 

lb. 

101. 

180 

131. 

19. 

185 

133. 

FORMAT 

41. 

190 

134. 

FORMAT 

100. 

195 

13S. 

FORMA t 

10b. 

iH 

136. 

FORMAT 

40. 

66. 

205 

137. 

FOHmaT 

114. 

210 

1  3fi. 

FORMAT 

117. 

215 

139. 

fORMAT 

123. 

220 

140. 

FORMAT 

12b. 

2<;s 

141. 

FORMAT 

lb. 

230 

14?. 

FORMAT 

20. 

?3S 

143. 

FORMA  1 

67. 

^4o 

144. 

FORMA  r 

bB. 

24S 

145. 

FORMAT 

98. 

66.   104.   119. 


49 

STATEMENT  NUMBERS 

DEFK.tU 

4 

MAXIMUM  NUMHER  OF 

RFFFRtNCES 

S9 

TOTAL    NUMHER  OF 

REFERENCES 

1.2 

AVERAGE  NljMHtR  OF 

RFFFRENCFS 
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ti    a«— oraaiajaa  — 
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uj       W  i-  z  zu 
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D  »-  uj 
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D  c  z  =)  «i  a  o- 
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i)  u   v  t>    t\   u  a 


-.  *  -O  —  ^43^-.^.^.^.^, 


i     .  i    a  ■  i » ■  i  r  i/l  r 


m  j)  ii  m  n  Jj   <£>< 


aM^,^^       ►-a— <aaaaarcraaaaaaaa 


«    -z 

1      3 


i/i 

i/i 

f 

UJ 

o 

o 

l_) 

z 

z 

z 

OJ 

.11 

a 

LI 

I 

UJ 

UJ 

" 

a 

at 

I 

U. 

w 

u 

<  J 

o 

<  ' 

L, 

■1 

a 

(I 

*? 

r 

-c 

I 

u 

J 

5 

3 

.iJ 

) 

3 

> 

^ 

2 

Z 

7 

■y> 

i 

i 

D 

t  i 

I 

* 

I 

5 

a 

O 

> 

t/1 

1 

►" 

0 

ir 

^ 

to 

IT 

•t 

ru 

Xi 

rv 

"' 
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d)      LISTING    OF     SUBROUTINE     GRID 


?4. 

?5. 


SUHROUTINf      GHllj      (Z,NH0wZ«NCOLZ> 

SUhrOUTine  c,mo 


PURPOSE  — 

TO  INlfnPOLAH  A  SfWItS  UF  IRREGULARLY  SRACFn  'hHEE-DIMENSIONaL 
DATA  Oh5fHVATIONS  TO  PROD^Ct  *  Ma'NIX  OF  EOUaIL*  SPACtC   (GRIODEDI 
FLECTIONS,  a  WEIGHTED  Mu«lNO  AVERAGE  IS  used. 


ThE  GCAhS  STSTtM 
PUM(,Ut  UNlVtMSlTY 


A.  kfITh  turnfRi 

CIVIL  ENGINEER!' 
JUNE  1Y6V. 


I1..AGE  - 

CALL  G"IO  (Z.nhowZ.nlOLZ) 
WHERE  -  Z  =  MATRIX  „E  GRIUUEU  VALUES  PRODUCED  BY  ThIS  SUnROUl  INF.. 

-  N«0#/  .  MAXIMUM   (DIMENSIONED  )   NUMHEH  OF  BOwS   IN  Z  MATRIX. 

-  NCOLZ  ■  maximum  (DIMENSIONED  )  Number  OF  COLUMNS  IN  Z  "ATHIX 
NOTt  - 

this  rout inf  utilizes  Labelled  common  /gdata/ 


DESCRIPTION  OF  CONTROL  CaRUS  - 

NONE  REQUIRED  »  ALL  RflRA-'E  r£kS  SPECIFIEIi  BY  CALLH-G  PROGRAM 


REQUIRED  SUBROUTINES  - 

SI  ATS 
OPTIONAL  SURROUTINtS  - 

NONE 


dimfnsiun  xiiooni.  Tuouu),  "(louoi,  ZinrowZ.ncolZ) 

DIMENSION    |){20>>    Nl<20>>     T  I  I  L  (  7 )  •    AI6>.    BU).    r  16) 

COMMON    /('.DATA/     N.Xt  T  .*,  I  1  Tl  •  ROWS  •  COt.  S«  lTYPFir,ST7E»HNUlLlS't|PUh.TA 
IPE'XmAXiXmIN.YmAX.TMIN 

ut»L    maxx.mINx.mAxy.mIntiMaxk.mINw.maxz.mIn^ 
INTEGER     lAPF>YtS>Bl>U.RO«StCULS 
INlEGE"     SYS!N.SYS0TtS0F«;<SllF3 
DAI  A    s Ys I N. S YSn T. SUE 2.SDE 3/5,6,2,3/ 
IE      <L I  ST .NF  . 1 .AND. LIST. Nfc .3)     GO     TO     10 
WRITE     (SYSOT.lZb) 
WHITE     tSYSflTiliO1     MIL 
DO    5     1=1 .N 
5  "RITe     (SYSOTil*0)     |iX(i>  iY( | ) iw( I> 

10  CALL     STATS      (N,  J  ,  1000.  l,A.XMAR,DUM.XSIG,MAXX.MINx) 

CALL     STATS     IN.  1  • 1000' 1 »Y« YHAH»nUM»YSIGi«AXY tMlNYl 

CALL     STATS      (N.  1  I  1000,  1  .W.KrtAR.nUM.wSIG.MAXW.MlNW) 

WHITE      (SYSOT.12?>) 

WRITE     (SYSnT,l«5)     TI1L.N 

IE      (HND.EO.ll     RO     TO     IS 

C         

C  BOUNDARIES    DEFINED    B'     UATA 

c  

XMAXaMAXX 

XMlMaMINX 

YMAX=MAX Y 

<«I';:M|  NY 

c        

C  BOUNDARIES     SPECIFIED 

c         

15  EN=MAXW-WINW 

WMAX=FN«MAXw 
WMIN=MIN«-FN 
MAX2eMAXw 
MIN7«MIN» 
M«0 


100 

1"2 
Id* 

I  (16 

ini 

I I 

112 

1  14 

1   1ft 

1  I  - 

120 
122 

1?4 

!?*> 

l?B 
130 
1  32 
1  3A 
I  »ft 
1  Ik 
UO 
U2 

144 
1  Aft 
!4H 
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40. 

4  1  . 


42. 
»3. 


4b. 

46. 


47. 

4H. 


so. 

51 . 


52. 

53. 
5*. 

55. 

56. 


59. 
61. 
62. 

63. 

64. 
6S. 


«.8. 
*9. 
70. 
71  . 
7?. 
71. 
74. 
75. 


76. 
77. 


FN«N»0 

DX'XMAX-XMIN 
DY»YMAX-YMI.N 
A«tA=OX»DY 
BESToSUHTfABEA/FN) 
AVEOIS»BEST»J. 073*6 

WHITE  (6»160>   XHAR|XSIGiMAXX,MINXiYBAM,YSlG»MAXY.MINY,wBAN,WSIG.MA 
),XW,MINW.AHEA.AVEUlS.bEST>GSlZE.RO"S.COLS 


C 

c 
c 
c 
c 

35 


c 
c 
c 
c 

*5 


C 

50 

55 


DEFINE  GwIO  SUE 


GRID  SPACING  SPECIFltD 


IF  (GSlZE.LF.O)  GO  TO  2U 

COLS=(1.5»DX/GSIZE) 

RO*S=(1.5»DY/GSlZE> 

M=l 

GO  TO  *U 


NUMBtR  OF  COLUMNS  SotCIFItO 

IF  (rOLS.LF.l)  GO  10  30 

GSI2F«U*/FLOAT(COLS-1> 

RO"S=(l.b*OY/GS!ZE> 

M.l 

GO  TO  40 


NUMBtR  OF  ROWS  SPECIFIED 
IF  <HO*S. IF. 1)  GO  10  35 

gsize=dy/float  (Ro»s-n 

CULS=>(1 .540X/GSI7E) 

M=l 

GO  TO  40 


NO  SPECIFICATIONS  GIVEN 
PROGRAM  COMPUTES  HESI  GRID  SUE 

GSIZF»mEST 

COLS=(l .S'DX/GSIZE) 

ROwSall .5»DY/GS1ZE) 

CONTINUE 

IF  (COLS.LE.NCOLZ)  GO  TO  45 

COLS»NCOLZ 

GO  TO  ?b 


WRITE  FINAL  CONTHOLS 

IF  (Hni).NF.l)  GO  To  50 
WHITF  (SYS0T.150) 
GO  TO  55 

WRITE  (SYSOT.1S5) 

WRITE  (SYSnT.i*,^) 

IF  (M.FU.l)  WRITE  (SYSOT.170) 

IF  (M.NE.l)  wRITF  (SYSOT.175) 

WRITE  (htiflo'  XMAX.XMINiYMAX.YMlNtGSlZEtROwS.COLS 

IF  (LIST.E0.2.nR.LlSl.FO«3)  WRITE  (6.1fl5) 

PERFORM  INTERPOLATION 

OY«(oY-lHOWS-l  .  l»Q!>IZE)/2.0 
OX"«OX-(COLS-l .)«Q5lZE>/2.0 

Ymax»Ymax-0Y 

XM1n»Xm1N»OX 

HGKID»GSIZe/25.0 

FN* (WMAX-MAXW) /2.0 

MAXX»MAXw.FN 

MINX«MINW-FN 

LOOP  Through  grid  Rr>w  hy  how 

DO  100  M»l,Rows 

YPwYMAX- (FLOAT (M-l) )»GSIZF 


B 

150 

R 

152 

R 

154 

8 

156 

B 

158 

8 

160 

H 

162 

R 

164 

R 

166 

e 

168 

8 

170 

B 

172 

& 

17* 

B 

176 

R 

178 

B 

Is  J 

B 

182 

B 

18, 

B 

186 

B 

18B 

6 

190 

R 

L92 

R 

19.. 

B 

196 

B 

198 

P 

20  0 

B 

2"2 

B 

204 

B 

?06 

H 

208 

B 

210 

B 

212 

H 

21» 

« 

216 

h 

2lH 

B 

22  0 

8 

??l 

B 

224 

8 

226 

8 

228 

B 

2^0 

H 

232 

H 

234 

a 

23h 

B 

?38 

8 

24  0 

H 

242 

B 

244 

8 

2.6 

8 

?4H 

B 

28  0 

B 

28^ 

B 

2S4 

8 

2S6 

H 

25rt 

B 

260 

H 

26<? 

B 

26. 

B 

266 

B 

268 

B 

270 

B 

272 

B 

2'4 

B 

276 

B 

2TH 

8 

280 

B 

282 

B 

284 

B 

2«6 

B 

288 

H 

290 

B 

292 

B 

294 

B 

296 

B 

298 

B 

300 

P 

302 

M 

304 
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LOOP  IMM'UUGH  H0«  COL  aY  CUL 


10  0, 

I  til. 


10?. 
1M. 
In*. 
10*. 

i  n6. 

107. 
Inn. 
In9. 
110. 

111. 
II?. 
lit. 
11*. 
114. 
116. 


117. 

11H. 


1  19. 

lao. 

i'i. 
122. 


1?3. 
1?». 
1?5. 
126. 
127. 
12H. 
129, 


no  <4b    J=l  .COLS 

XF=XMIN»(FLO«1 (J-l) ) •bSI/E 
OU  60  1=1,9 

60         r>U)=i.Eio 

n1 (ll =1 

Nl (21=1 

1)0  HO  1  =  1  ,N 

nisi  =  i  (xp-a  1 1 1 )  "2)  •  i  (yp-v  1 1 1 1  ••?) 

IF   (OlST.Gt.HGKlU)   GO  IU  6S 

c   — - — 

C    AN  OhigJnal  OAlA  POINT  LTFS  WITHIN  l/?5Tn  G«T0 

C    ShaCInG  .THfcPFF  oPf  INTtHPnLAT  IUn  sToPPfco  FOR  TM|S  POInI 
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??h. 


2?7. 

??H. 
??<». 

?3o. 
231. 
?i?. 
?33. 
23*. 
?35. 
?3*. 
237. 

?3fl. 
23<J. 
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^\ 

s  r*  r 

I-.H1I          HI      M 

Hi 

Utl    I.-Jt 

J 

s 

I..AA  1 

or 

KL(  h  Kt 

ict  s 

db 

I'JI  .11 

,.Mm-« 

!!(■ 

Htl-tHt 

<cts 

1  .i 

-VM- 

I  > 

ntfHt- 

'Ct  s 
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I")  f» 


f.    AI  >P*-^**^J 


-4       <v.  a  a.  rv  a  — 


»        aj  aj 


r\^-Tif  rcc    i 


•a   t  u  u  id  i 


-    i         _',_._j3r        —    ,  — 

l  — •  Aj    -J_     —       >-  —      5        5. 

i  —  —  y  7   t  t   x  y  ir.  in      »-  nj  *, 


tit 


«J  I  L 
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j)    LISTING    OF    SUBROUTINE    STATS 


7. 

8. 
*  . 

in. 
11  . 
l?. 
13. 
1". 
15. 
lft. 
]->. 
IP. 
19. 


SUHHOUTINF  STATS  (HOkStCOLS. I  I l JJ.X.XBAft.XVAR.XSIG.XMAX.XMIN) 
C.  SUBROUTINE  STATS 

c  . -. ... .... — ....... ... 

C  PURPOSE  - 

c  ro  colcui  Ait  statisiical  measures  of  one  or  t*o  dimensional  data 

C  ARRAYS.  "tASURES  DETERMINED  ARE  1 ) Mt AN , Z ) V AR I ANCE . T ) ST ANDARD  DEVIA. 

C  T  ION. 4  )  MAXIMUM  VALUt .S) MINIMUM  VALUE, 


USf&f-  - 

CALI 


STAIS IROhS.COLS. II.JJ.X.XbAR.XVAR. XSIG.XMAXi XM1NJ 


ROWS  -  NUMHFR  OF  HOwS  IN  ARRAY  BEING  ANAL  YZEO .  (  I  NTt  GER  VALJEI 

CULS  -  NUMf-F  R  OK  COLUMNS  IN  ARRAY  BEING  AnAL  YZED .  (  I NTE  GER  VALUE) 

II  .  MAXIMUM  NUMHtR  OF  RO*S  ALLOntO  FOR  ARRAY  BEING  ANALYZED. 

jj  .  maximum  NUmhER  OF  COLUMNS  ALLOWED  FOR  ARRAY  BtING  analyzed. 

X  -  DOIA  ARRAY  HE.  INh  ANALYZED. 

XhAw  -  MtAN  OF  ARRAY. 

XVAR  .  VARlANCt  OF  ARRAY, 

XSIG  -  STANDARD  DEVIATION  OF  ARRAY. 

X">AX  -  MAXIMUM  VALUF  OF  ARRAY. 

X*IN  -  MINIMUM  VALUF  OF   ARRAY. 


DESCRIPTION  OF  CONTROL  CARUS 
-NONE  REUU1RFU 


REQUIRED  Siihpout  InF  S-- 

-NUf'F 
OPTIONAL  SUBROUTINES-- 

-NO«F 


INTEG 
ni«EN 
FNO«F 
xS'iRi 

XSuMl 
XMflXl 
XMJNi 


ER    ROkS.rOLS 

STON     XI II ■ JJ) 

nwS«COI  S 

n.u 
n.o 

X(lil) 

xiI'D 


nn   s    1=1 .rows 
no   s   j» l. cols 

IF  (X  I  I  , J)  .CT.XMAX )  XMAX"X(I,J) 
IF  (X(I ,J)  .LT.XMlN)  XMIN»X(I,J) 
XSUM«XSUM.X ( 1 , Jl 

XSOPaXSfJR*X  (I.J)»X(I,J) 

XHi.RoXSUM/FNO 

XVAR*(XSOR-?,n»XHAR<»XSUM»FNO«>XHAR»XFlAR)/FNO 

XSlGrSORT (XVAR) 

RF  TURN 

END 


STA 

2 

STA 

0 

STA 

6 

STA 

a 

STA 

10 

STA 

12 

STA 

1» 

STA 

16 

STA 

18 

STA 

2n 

STA 

?? 

STA 

?<. 

STA 

26 

STA 

?s 

STA 

30 

STA 

3? 

STA 

3* 

STA 

36 

STA 

38 

STA 

»0 

STA 

»2 

STA 

M 

STA 

4b 

STA 

»8 

SIA 

SO 

STA 

s? 

SIA 

54 

STA 

Sb 

STA 

s« 

STA 

bU 

STA 

b? 

STA 

b4 

STA 

bb 

STA 

68 

STA 

70 

STA 

72 

STA 

7% 

STA 

7b 

STA 

7H 

STA 

80 

STA 

82 

STA 

84 

STA 

8b 

STA 

88 

STA 

90 

STA 

92 

STA 

94 

STA 

96 

STA 

98 

STA 

100 

STA 

102 

STA 

1  0" 

STA 

10* 

STA 

I0B 

STA 

110 
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k)   FORSTAT    CROSS-REFERENCE    LISTINGS    FOR    SUBROUTINE    STATS 


SlMWdll  I     t     s  I  ..  I  x 

si  An  Mf  -,|    MuMHt  u        at  ^  if'M  Ion 


RtFtwtNCFS 


1     SrMF-tNl    NUKbEWS 

OFF  INK) 

?     HjMXU"     NU"HtH     OF 

MtFFKtNftS 

J>     Ti'T«L           NUMHtH    Ot 

HEFtHtt-'Ctb 

2.«    a\.f><A'-F    r.n-'i't  h    nt 

I'FF  1  "i  'ft  ^ 

SUMRoUT  pif       Sl»  IS 


LAHFI 

trPF 

OFF 1NITION 

HEFtwtNCES 

COLS 

1 

2. 

4. 

10. 

F  NO 
I 

u 

1 

lb. 

1 1  . 

1ft. 

12. 

1  3. 

1*. 

II 

I 

n 

j. 

J 

1 

X 

1  )  . 

12. 

1  3. 

14. 

JJ 

I 

D 

j. 

HOWS 

! 

<: . 

*  , 

9, 

SQHT 

p 

17 

.                 FUNCTION 

X 

u 

.       D 

J. 

7, 

B, 

11.          12 

XAAR 

P 

lb. 

lb. 

XwAX 

R 

1  . 

1  1  . 

XMIN 

U 

b. 

12. 

xsir, 

R 

1  1. 

XSCjR 

R 

1*. 

lb. 

XSUM 

R 

13. 

15. 

16. 

XVAH 

R 

16. 

1  7. 

lb 

Sv^mOLS 

OFF INtO 

7 

MAX [MUM 

WIM.ItH    OF 

REFERENCES 

19 

TfTAI 

NUMI1FH    OF 

PEFFKENCES 

2.4 

Avt  ..fll.F 

Ml  <"tf    OF 

PFFFMt  NCtS 
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I)     LISTING    OF    PROGRAM    CONTUR 


I  Ltn  I  I    .fUl   i  OOlfUl  ,  t  I     .1  <  I  /'Ft    .1,  T  Af .  '.sINhJI   .  I  Af  F 


f ,  ul     rui.l 


.UH     -inks    Cf     l.h||,UfcLi     I'M  n     [li\ 
ibtU     11"     A     M   H  It   ,     J      1[      PH|I,.K«' 

•  ICnli.A'     I  F  Ha..  I  '-'r  i.  I    OF     (,r 


'£      CAl  Cuf  f     M|. (lilt' 

•  n  1  T  I  f  i.    ri  I    f  -  a'-* 
-..Pur , 


InF     uC. 
futojot 


<b    brb 1 F  m 
I.JlVFf.1  1  I 


A,  M-  I  In  lu.M-  h 
(  !  v  U  t-Mull  Ftnll 
Ji     »t      l-»h». 


ii("tM»l      I'. bit 


-Mb 

InF 

PL 

ru 

m     » 

ll .... 

<p 

.11- 

1»  1 

■  II' 


1AM     Cxf.     tVi'U"LF       MAWS      111      Afj  Y      SC«LF       KllHlM     1ft     LIHIIS 
"1HH      «|L/1M      (I    A*1"UM     e      ,-H      lt/C"FSI. 


tt  Ll'«     Ir 

■tttCIS. 


ollliinb    ant     Ai/Al|i.-Lt     r,  4(J    Ci'i    Mt     s^triHh      lh«'u-n     Iff     us' 
.      l.Mllb.      Int       ,SFH     ,..,   t     b.tCltl      (0      ll)-,i-S     Al      HhbULl,n     Uh      IxHt..- 
.  Vlil    SltNII      -IAf       I  I    A,'bF  LIFim      ■OIF-     (.OOHI'1  i...  I  F  b     hi       InA'.SL,.  IlUN     liH 
,nl,ul     "Nlr     ilF      I  rlt      I  nt  f  t         AJUH     4HFS     I     >,li(m     /I. 

'<U."«»     CA.     f  rt(i|,U(  t      IbdrtlHK       Hl_"0      UlAbHAt'b     OF       lilt      CO  '  I  U'Jt' 
ilMllMl)     Int     UAlA     (,m.u!      .  LL      I  '•  -<l  F      »ttS    utMInt     PLIllTINu.      Int 
-ti-->f  t  CI  I  »t      Jfl]Ui       )N    (     ■  "i «  1  i .  i.  I  |n        .  11/1     lkAl.bt.Al  IU'i     -ILL 
r-L  ft  I  I  I  imi,    Ul-     IFFl     AMI    hll.nl      v  1  (•  wS    f,,h     b  I  I- Kt  oaCOf  If      I      i-lil 


I ->f    -kik.kam    )b    Of/nLF    il-    t-i'iillf.--    a    Nu-.ntn    01-    nllimt.il     i-Fb 

'IF      »     SI-lftLt      '■Fl      OF      LIAl/li     ti'Cn      Ihi/DLV  |M.     PIFFtnt-.l      M»F     b(  ,  L  t  S  i  C'">   '  I  0' 
l(     IF  KVULil  Ml' I  "  '  lUlb     U<      It-A      S(   A  I  |     tub,      «llfv,U|      Ft(ntAl)|Mo      l-t      l)AIA. 
Kl.LUtl    OF     CLI     llJUrt    LIM:     L00m|J)>ia1f  b    Ant     CaiC.ii.aIHJ     A-mJ     b  I  f.MF  t)    UN     A 
V"»lr-i    I    f  t  .     InfcSf     vuf-'tb    „t-t    imuIMHi    a    U    fLO[T|-u    F  uh    u(n    m„p 
■*FU')tblni    bO     Intl     IrtF     CAl  C"L/  I  IIj.iS    A*F     Ml     I-WFu    A»-u     Iff     PLOT  I  I  mu 

'IF       F'OL  I   I  ft  t  'FS      lb     tFF    1(    IF  1,1  . 


bL'.ilCn    )«Pt     I     l'-fl3    I     lb    obtli    10    blimt    COi\  I  OuH    LINt    COUhI 
-l    •   IF  b    F-HIOH     |i,    FLI-1  Tl»-G. 
f  L  ■    I      1  Aft     I  3    ff  liUll-fclJi 


OF  b(  f|f  I  ]  'I. 

SfF  C   I   ' 


F     CO  .  I KOL     CAHOb    — - 
.OIF     - 

•'/L      (HMO     CO'.l/.lr.b     A     CONMOL      ltlFf.ff      ("CON)       IN     C"L"fNb      !-£•, 
II       .b     KtlJ'Jint       plfUl      IIAlA      1*      IMF      F  0W1     >,F      V  ,.  h  1  A  n  I   t  V     ANO/I'H 

lulUf'oAL     LA'-'ub.     (iffi      X/.'-UhLtS    MfUUIKfU    bni,UL(,     ..fc     PlACtl 
■>      I.IMI-«,^.)^1J.*/,     tlL     «|Irt    -FriMAL    fOIMb    PUNC'ifU, 
Ul    I-miL     lulF'.F"     l'[ilC»IIS     Ihsl     «r<Ul  I  I'I'iAL     (.nnl-S     n.-t     Mt  CF  bS- 
(.(     fltlF       'f  F'<i,l  lnl,,      Intbt     CA"Lib     briOOl  II    |I)ILU»     t  M»  t  L>  I  A  I  f  LY  • 

n F  MAFln  S 


ll       ,-    MFAU     A  4|,    C"1"'1"1^     oaIa     "OlfI« 
ML  '  b-  1  I  NLMhfcH    Of     f|   »S 

i)    iO-'Mt  '<    OF     Cl  LOM  .b 

3)  .  [Din    i,F      -"F     Ii»     1  ,(  nts 

<tlbC»Lt     FnCIU,.    In.     /     VAI.UFS     l/«///SColtl      ilF     KL«''«»I 
Ci'Ub-     Int      J  a  I  n    nl-I»     lb    ftou     IN    '«l         BV    mj«    Fmj.     Igu     lu 

n'HToM     I     FACm    nut      IS     bfANNFfl    F  h  ,..     LFtl      to    f|i.fl). 

I  Fit       -'AIk|.      |i    UFA!'     AtClfl'l'l'     10     '•     F-ft  «  I  I"  iH  t     Itl-l'-t' 

Fi ■>  I    istt     ,r  iii.=-i  si  I  .    't>i    ■»■■    si/t    i^    !(.->. i'"    . 
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IIJfH     I  {■<  --Klllfl  It      n>f      v  1  F-  •  I  r.. .    Plblll'l        '.IF     1"E     Su-FA(F      I  "t  I  1 
DE'-Hcts    !'■     |r^     *      ■MtCl  low. 

/A^I-.rlLF  5-1  )      InEU      IN     UFfjHEFS. 

Hi)i|L     CA-ub-       NUNC 

IPFFI,.  U0 m-KOTAIE      ImE     '/IFbINn    PoSUlO,    OF      1ME     SuMFAfF     Intra 

DtCiHtfj     IN     IMF      *     NlHtC  I  ION. 
ja,  I  .,(3LF  b-l  )      Intlu     [N    UFOfctS, 
noOL    O-mS-      nvJNt 

■  iHl-  HA  [  If,  ...mlU  It     Inf      >U«II«h    PUSIIIOp.    OF      I  nt     SuKrACF      IntlA 

LiENrfEEb    IN    InF     ^    nlhtC  I  ION. 
«Ahl,lf)Lh  b-1)     IrtEU    In    UFuKEEb. 
AUDL    CA"i)S-      NOnE 

DPFK.,  I  Ioim-CALL    SJiiKciUllft    f-JUJ    FO-<    USEP-litFINEL    '"On  I  F  I  CA  I  IUNb  . 
i/A*I."nLFS-    NUNt 

IDOL     C«-l'S-        NUNt 

DPEK.l  I  U'll-PLOl      03EHS     flift     AWL)     ACCuUNT     Mlr-.HtM     ON     PLOT. 

JA-.I   mLFb-    NONE 

(ll)OL     CAmUS-     NAi«E      AM,     ACCOUNT     NUI»btH 

DPFXjIK   »-IDEnIII-V    CO'.Icuh    LtvtLb    Ml    PLOTTING     I  Ht  I U    tLEVAUONS. 

KAHlArtLFb-1)  HE  It>ll     OF     MMHt*S     UN     iNfrits     I 

AI1DL    Cc-i'i-    «     iu-C"AH6CltK    F'JHmaI     In    CUI.S  1  -  I  0  .  g  1  v  1  hr,    fuf-maI 

FOH    r.Lf  VAl  KuMb-    UbtU    'o    sPtrlFf    NUnptH    of    IUCMaLS. 

DPFrt  i  I  l0,«-PLUl     COnIOUkS     »  I  1  m    PtKbPtCUvF      TWANSFOKI-lAT  lUN     10    PNOt 

-OCt     LtFI     0-     kIGh(     Sltl'FO    -OiJFLS. 
■/.WI.idLF  S-l  I  I'hlmCt     Ml    PtniHtC  Ilvt     PLA'it. 

?)l)IblA'*CI:      10    OhJi-CI     PLANE, 
5111)1,     CA  -US-         IMJ.Nt 

dpfiVuI  liiN-lrtANiLA  I  Ion   of    co.,iuum   cuOwOIna  !(■&. 

VowIahLFS-1) TWaNsLATIOn  OF  «  In  InCHFSIiOT  SC  AL  E  0  i  I  !■  UF  KtAauHF) 
?) iHANbLATlON  OF  y  In  InCmfSIiOT  bCALt  u.  1  F.UF  rfbuHf) 
3IIKANSLAI  Iun    OF     /     In    InChF  S I   .oT    SCALEUiTPUF    htAbijHE) 

ftDOL     CAmlb-         NUNt 

llffHA  I  Kill-     OtFi'»t      KlOllFH      luKt      Ab     0  I  Ht  K      ImaN     NUKP1AL     NUMhEHa 
iMHl.riLFb-l)      PLOITfF.      TAFt      NIMhtH 
A')i)L     CA'-llb-         NUNt 

opkh.iII'.j-    ritAu     t  1  1 1  t     K'H    Pi  ul  (PLul  TFO    J     TOP    uF     h*P|, 
^AnInHLKb-DDIbrrtNCf     FHfiM    LtFI     tOlit     OF     Kt  fi  I     TO    blAPT     (IF      UlLt 
I      IN     [NCMf b     )  . 
•  iMElGnl     UF     I  I  I Lt     Lfc  I  lt*b. 
AOOL     CAwllb-     IIILt     PUNCrltU     In    COLUMNS     I  -  HO 

DPFrM  HON-     TFHM'MAIf      PL<    I      T  UPE 
VAW  [.,dLF  b-     NONE 
UDDL     CA»Ub-         NUNt 

OPFWaIIiin-    UtAj    LAnFLS     AND    COOHU 1 NA IF S ( "A > ■     PF*     PLOT     3     Jul. 
^AM|AhLFb-ll      *     C0UW|llNA1E     OF     LAnfcL 

?.)      1     C'JUHl'll.Alf      l,F      LAlEL 

i)      I     C'JUhIiIi    Alt     OF      LAntL 

-.)      Mtl'jHl      OF      I   AUtL     Lt  I  IFHS. 
SOUL    CA'-ub-     A     LAltL     PlJNCMtL/     IN     COLUH'S     l-T«i 

i)PF.w.i  I  In-i-OEF  Iiit    KH.UiAfLY    SPACEO    COjIOiihS. 
VawIahLFS-1)  bfBH  I  Ini,    COi«TOLiP 

?-\       CON  rOUP       [NlFHVAL 
J)      F   InAL       IH  |c,hE  b  I  I      CUlvl  OUP 
AllilL      CA^OS-         NUNt 

DPFrtA  I  ir.N-kEHEFlNE    COnTOUM    [NltPVALS    FOw    Samf    uaIa. 
VAMIi'UFS-     SAME     J,    F  uw     'Cur^lb 

i)F>F  U,<   I  K  n-UEF   lot       lF>Mt';U|    AHL  T      SPaCFU     CONlO.JWb. 

tf»wT  ..iiLF  S-l  )  uCun  I»mwmEf    OF    COnIO'IHS. 

Al)l)L     CA'ijb-1  I  F  OH  ■  I  A  I      OF      CONlmjHb 

PICoNMiup    LtvFLb    AS    DtFlNF.O    am"VE. 

'WiTAlli.-w     SF-fcCIFC     ninth     OF      PAHEHI      11      (]J     ?H      iNCHtS) 
Jli-ll  ..fiLF  b-1  I      *  I  0  I  n    nF      FjfK, 
AHUL    CA'-oS-       MONK. 


I  h  I 

1  h? 


1  1? 


I  JO 


?1  ? 


??H 

?3n 
?1? 


?il 
?4A 


?b>* 


?F|A 


a 

?H? 

A 

?H4 

a 

?««, 

A 

3  J  H 

1 

?<<ri 

a 

v  )? 

A 

?<*4 

A 

P  J« 

a 

J  <- 

a 

inn 

A 

T|? 

a 

^1  "A 
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I.  >.'►-(.  11,    .-    SffcClrf    MJwm/.I    Of    OmIA 

*..w|       ll  V  S-     *l)«c 

1III1L    (."•  Ob-F  "h,i»  I    of    r,oiA    maThIa    I  N    tol   ,•  '  S    1-hO 

?,  OPF*  J  I  liw-PLOl    »jM    nCCOi  UlNii     10    ALL    COnThi>|.S    GIvcn    M'-Cf    laM 

CALL     [J    Hl.Ol     ILfcbl    NCON«20) 

J  w  |  „-iLl  b-     NONE 
AOOL    C»«  JS-    NO  «e 

PI  JPFP.i  111,  <-     C"A  "lit      »Il)lM     UK      MLUIItl)     "HP      (CHANGE      bLALM 

i/ArfliirlLf  S-l)     NEW    »IDlH    IN    I'CHES 
Al)l>L    CABUS"    iiOimE 

CommE  *Tb    nN[l    HtStHlCI  IONS 

PH|Uh  Id  HEAOInG  I  >N  Uala  MATP1A.  IhR  fn-««IINCUN51VHi',li  COvIUU*- 
II-.rFPVALS(NCrtN=15    Orf    17)     MUSI    PE    DEFINED. 

tFFo-E  ruLL  [0  PLUI  (ImC0n=2ii)  FurFP  wIOTH  INC0N=1U)  A*u  ACCfU'il 
CiiDE     (NCiKsuM     iusr    HL    OEFInFO.    PLOIUnG    1hF    r.  it  h    PEbtlb    oLL    COmIi-ulS 

[')    PiOl     u    Nt*    it  I    UF    data    USE    NCliN*l     ill"    OH    *lInouT    a    Nt« 

Fiih-aT      UK     CU      fOuH|Ni,     COl.TMOLS. 


I'Eb 

THE    FULL"«IN<,    ENTPY    POInTS    - 
NUMHfcP.F  «(.  rUK.WHEPF.  ,L  tl"t 

ni>a«D    CALrOMK    HOUT  TIMES    6b    DEbCxJufu    IN     I 
tPT    Few    NOHrlEP.    «imIcm    "A1;    bEEN    MOOIFIFU    A 
LOWS    — 
A.Y.rtF IGnf.NUM, ANGLE .  FOPMAT ) 

sTaHTIng  coukuInates  of  The  Lu»t*  le f  f - 

NFH    of     Tut    FIHST    chaHaCTEH    In    i liF     STPIno 
llESIPEl)    mEIohI     In    INCHES    OF     Tnfc    i^hapaC  ' 

UNSlht.Ti     JAPIabLE     IDE.'IFIEP.    up    FAPHfc.bb 

bf     vAiOE     lb     lu    ,1E    PL0T1EU. 
ANGLF     1  (^    UEukeES    rOUo If nCLOC««IbF    FMOM 

DIHFCIION     Al      WHICH     ThF      KUMHEPS     «KP     70     tit 

ITEO. 

OLLEPIIh  STWlNG.  ni  HOPE  Than  io  rnAPAC  1 
luU  IHE  FOPTKAIn  IV  KUHMAI  Hf  »H1Ch  NUM  I 
CONvEPTEU  FOh  PLUTTIN,,  I  USEU  10  RIVE  NO 
IJFC  lr',/>L  hi  ACtb  )  . 


10 

I' 

ti   ■ 

HP 

Si 

11*00  r 

I' 

CI 

u,  1 1  wr, 

1  s 

Si 

MrlOL, 

bt 

a> 

E     bU 

OALS 

.    E*C 

DUE 

AS    F'JL 

L    ' 

iUl 

Ntrt       ( 

nO 

r 

-    The 
CUP 

,M 

-     (rlf. 

-AC 

XrtO 

le 

PLU 

t»b. 
ION 


Of 


PUPI1HE     M«i     IMPLE'iNTtU    rt    COPYPLI     Llnw^Pf     hOOlINt     *HlrH    PhEPAMEb 
PLOT     [^tS    FOK    PLOITInG, 
■I)     UIHt"    KUUTIwES 

-  SC  i  i    -bCANb    OAlA    TU    Lt'CAlt    CUNIOUPS. 
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10b 

104. 

i  15 

1  1    7 

1  01 . 

12" 

i  n  m 

1   "t;  . 

12S 

l  i  i 

21. 

13' 

i )  ' 

21  . 

i  Jb 

I?  i 

21  . 

!<,,! 

l?9 

21  . 

14S 

112 

21  . 

ISO 

1  *b 

21  . 

15b 

U  1 

13b. 

lhO 

ISi 

1-/. 

I6S 

IbJ 

lbu. 

1  I 

l"-h 

lbj. 

i  rb 

15S 

1MM. 

180 

1*2 

142. 

1Mb 

IftH 

167. 

190 

1  74 

173. 

1  <Jb 

17b 

173. 

20C 

i  rs 

167. 

20b 

IRQ 

17b. 

210 

19b 

21. 

21b 

19H 

13b, 

220 

<;rn 

H  . 

225 

202 

4b. 

1  it) 

20« 

162. 

23b 

2n6 

20. 

24' 

2    7 

199. 

2*5 

21  9 

,                 FORMAT 

lb. 

250 

21  0 

,                 FORMAT 

19. 

25b 

21  1 

,                 FORMAT 

7b. 

2*i() 

212 

■                 FORMAT 

72. 

265 

21  J 

•                 FORMAT 

Hb. 

27- 

21* 

.                    FOR 

-AT 

losi. 

'7b 

2lb 

.                    FOR 

-AT 

19H. 

?H" 

21b 

i                    FORMAT 

200. 

'Mb 

2]  r 

,                 FORMAT 

20,2. 

g9n 

JIM 

FOR 

<A  | 

204, 

?95 

219 

.                    FORMAT 

20b. 

30" 

2?o 

i                    FORMAT 

207. 

REF  tHl^Ct-S 

6H.  7i).  73.  77.  h4.  H9.  9i. 

11U.       lib.       1?B.       131.       13*.       19b.       197. 


201.       ?o3.       20b. 


9S.       126.       132. 
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«J 

101AI            NUmrFH    OF 

REFERENCES 

l  .» 

AVt-Wit     MIIMHfR     nF 

REFEmEnCFS 
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n)      LISTING    OF   SUBROUTINE    SCAN 

Sum.  )  i  r  |«<t     ->f.  AN  •< 

C  H 

C  bom  uuIINt    bCau                                                                               H 

C         - - --       i 

C  H 

C  PllHHUSE     —                                                                                                                                                                                          1 

C  ro    SCAN     rut    L)AIA     AND    lOCATt     Tut     CONT'ILMS.                                                                            r 

C  X 

C  — - H 

C  THE    uCAHS    bvbTtM                                                                                                        R 

C  PUWOUE     UnIVEMSJTY                                                                                                                                  R 

C  " 

C  A.    KEITm    Tu-NKh.                         t* 

C  CIVIL    ENoInfEkIMj.                    R 

C  J"  *E     lVbv    .                                       R 

C  — --        H 

C  H 

C  'ISA6E    --                                                                                                                                                                     H 

C  CALL     bCAi'l                                                                                                                                                                              H 

C  H 

C         — H 

C  H 

C  neSCHIPTIun  OF  COnTHOL  CAHQb  -                                           R 

C  NONE  HfcQUIHEU  -  ALL  NECFSSAhY  KAhAME TE^S  PAbEEU  VIA  BLANK  COMM    R 

C  H 

C    - H 

C  M 

C  HEQUlHEo    bUhwoullNES    -                                                                                                                       R 

C  -    TrtACt                                                                                                                                                           R 

C  -    CAlC                                                                                                                                                             H 

C  H 

C  OMION4L    SUHHOUIInES    --                                                                                                                             H 

C  -    NOivE                                                                                                                                                             m 

C  <j 

C  — R 

DIMENSION  AM(lon.lUU).  KEC(bOO).  »lh00li  Y(lbOO).  lPI(3iili  1NX(H  R 
I).  lNr(M)f  IOH(?0li  AL(lO).  ailOli  ZLI10).  1  Ad  (  l  o  t  b )  .  UAI2).  OT(?  H 
/),    l)/(2)t    kTUOIi    CFhETAC),     sFmEIAk.),    ZUSftOli     T I  I  f  u  i  .    «OHr<  (  1  0?*       R 


9. 
10. 
11. 
12. 
1  J. 
14. 
lb. 
lb. 
17. 
18. 
19. 
21  . 
21  . 
22. 
23. 
2*. 
2b. 
26. 
27. 
28. 
29. 
30. 
31. 
32. 
33. 
3*. 
3b. 
3b. 
37. 
3b. 


ai 

COMMn 

I 1NX.I 
?MAA,C 

UX=1 

NP=0 
OL=D 
N  T  ■  M 
m|=HN 
IF  II 
IPFI1 
Iff  (1 

icrn 
i  p  r  i? 

IPF<? 
IPT  (  3 

IMF  (1 
IPI  (3 
lNA<  1 
INA(? 
[NX  (3 
INA  Ct 
INA  (S 
INA  (b 
INA  (  7 
INA  (M 
INY  (  | 
INI  I  ? 
INY  I  3 
INY  I* 
INY(b 
INYIb 
INYI7 
INY  (H 

I/x«n 

A  I  w*  F 
DO  lb 
MECIJ 
I5b«n 


FiNT.NI . I  A, IY. IUA. IUY. ISS. IF . IV.NP.U. jT i PY . ME C , CV . A i Y. IPT . 

NY.UL.AM.I0P.0P.DOlILAd.AL»YL.ZL»L4a>0A,0Y.0/,AMA)(fMT|SC«LE.Y 
Trie   IA,blM£7A.FM3,FNNlNN.M,CLtU»NLlN£b,NCUHV.Zl I IT.kohk 


U)     in.  10. b 

.  1  I  atl 

,21  =  1 

«3)=i; 

.11-7 

.31=3 

.  1  I  «b 

•  ?)»b 

.3)  =« 

)■-! 

I  =-1 

)  » 


)«-l 
)"-l 


I  OS 

in 

i  ? 

14 
|A 

I" 
I  21 
122 
i  ?« 

I  ?b 

I  j>H 

130 

3? 

3« 

it 

jr 

14(1 

1*2 
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39. 

40. 
41. 

"i. 
43. 

**, 

45. 
46. 
47. 
48, 
49. 
50. 
51. 
52. 
53. 
54. 
55, 
56. 
57. 
56, 
SSI. 
60. 
61. 
62. 
63. 
64. 
6b. 
6b. 
6?. 
6b. 
69. 
70. 
71. 
^^. 
73. 
74. 
7b. 
76. 
77. 

7H. 

79. 
HU. 
HI  . 
B2. 

H3. 
H». 
65. 
Hb. 
HI. 
8H. 
H9. 

90. 

'i  . 
92. 

93. 

94. 

95. 
9b. 
97. 

9H. 


*i  i«ir-i 

>U    ill    1  =  1.  -Tl 


Ob 

70 


jn    !=1.-T1 

IF    (4M1I.D-CV)    20t30f30 
II-     i^mi  1.1  »1>-CV)     30.25,25 
I  *  =  I  »  1 


I"  =  l 

IUX«-1 

inr=  i 

COLL    ThaCE 
COnTInOF 
f.T  1  =  f vi  ?  —  1 
uo    <,s     1  =  1  .Nil 


41    I  =  ltivTl 

IF     ( AMlMftll-Ci/)     35,45,45 
IF      (MM(f f • I » 1 )-CV)     45,40,40 
I  «M  I 


ir»l*i 

I  OA=U 
I  , )  T  =  -  1 

Call    inset 
CONflNiJF 
llU    60     1  =  1  ,~T  1 

Mf?""1 *•-' 

IF     (, 

IF 

1 

I  .- 

Ill*=  | 

Iiiy  =  n 

CALL    THACE 
CUNT [NUf 
l)U     75     1  =  1  .  Mil 

!MT2=mT*1-I 
IF     (--'■-■- 


ir?»ii  «i-i 

[F     (.i,4(»l2fNI  1-CV)     SO, 60, 60 
;F     <am^T2-1,M)-CV)     60,55,55 
:<=MI2-1 

It  =h  i 


;mT?=mT»1-I 

IF     («M( I ,Nl2)-CV)     65»75,75 

IF     (AMI  I inT2-1)-CVI    75,70,70 

IX*1 


I  f =NI2-1 

lOASil 

I I  >  r  =  i 

CALL     IhiCE 
CONl iNUt 
155=1 
N I l=Nl- | 

mt  i<mi- 1 

no    los    js2.NI 1 

uo    i 


los    .1=2, Ml 
105    I  *  1 .  *i  r  i 

IF     (/IM  (  I  ,  Jl-CV)     6(1,105,105 
IF     (am(  1*1  •  Jl-CV)     lo5,b-5,H5 
CDM«|00»( 1*11  • J 
IF     (  .HI     VO l 100 .90 
110    9S     lils  1  ,NP 

II-     IwtC(IOl-COM)    95.105,9  = 
CON  I I^UF 
I *  =  1 »  1 


IV=J 

in*»-l 
1 1  >  r  = 


CALL  fkflCE 
COM INUF 


R 

20  n 

B 

?o? 

B 

?n« 

H 

2  0  6 

R 

?OH 

B 

21" 

fl 

212 

H 

214 

8 

216 

H 

?  IB 

fl 

220 

H 

222 

H 

224 

H 

226 

B 

22H 

fl 

231) 

B 

232 

H 

2  34 

R 

?  ifi 

R 

?  |H 

H 

?4fl 

R 

242 

ri 

244 

Ht  TUHN 
ENU 


H 

252 

H 

254 

H 

256 

R 

->SH 

•1 

260 

H 

26  J 

H 

^4 

B 

266 

i 

^6B 

B 

270 

B 

J  f  J 
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LISTING    OF   SUBROUTINE    TRACE 

c      c 

C  5UHHJUT1NE     FHACE                                                                              C 

c         C 

c  c 

C  Pl'hFOSF     --                                                                                                                                                                                         C 

C  10     ThhCF     OU  I     EmCA    CONTOliH    LINE     IW     rtlMN.                                                                                  C 

C  C 

c      c 

C  Inf    iiC«Kb    SfSTEH                                                                                                        C 

C  PUHUUE    UivJVEKSIIY                                                                                                     C 

C  C 

C  A.     KEI IH    TUKNtHi                             C 

C  C.IvIL     ENOlNEEHlNli.                         C 

C  J'INfc      1969     .                                                 C 

c      c 

C  C 

C  us»GE    --                                                                                                                                                                    C 

C  C«LL    TNaCE                                                                                                                                                  C 

C  C 

c      - - c 

c  c 

C  nt-SCHIPTlUN    Of    CONTriOL    CAKQS    -                                                                                                           c 

C  NOME    Ht'jUlHED    -ALL    nECESSaky    PAHAM£t>rtS    haSSEO    VIA    BLANK    COMMOli      C 

c 
c 

c 
c 

L 
C 

c 
c 
c 
c 


rfMJUIWEl)     SUtlWOOtlNES 
-    CalC 


Optional    suhhouTinES    -- 
-    NO.lE 


IMMERSION  ami loOilOO) •  ki-CldOO).  A<1500)»  YI1500).  lPT(3.3)i  INAIh 
l)<  |Nri«).  IOP(?oi>  KL(lO),  YL(lO).  ZLMOI,  LAd(lutb).  UAI2).  0T(? 
•  )•  0/(2li  hTIIUIi  ClHETA<4).  MHtlAUli  ZI|5oO>.  IIIIBIi  *OHA I ] 0?4 
3) 

common   ii.^T.M,M,iT,iUAiiuY,is,5.n,iv,NH,i,,ji,pr,MEc.cv.».*iiPr, 

1  INA.lNY.OL.uM,  IHHtUHiUOt  ILAH1A1.  .  YL » ZL t L Ad . OA , 0» . DZ . AhA < . HT . 5CaLE  .  Y 
^MAA.CTHr  I  A,srnETA)Fii?.M»*''/tr'N.M,CLiU»NLlMES,i.CoMV,Z.  I  I  I  ,  WOrtK 

COMMON     ."AAX.MAAY.  11  N  A  1  Wl  NY,  SIZE 

HEAL    mAaAimaXYimInAimInY 
PY=0.0 


7. 

JT»II 

a. 

N«0 

9. 

IA0«I» 

10. 

110=11 

11. 

ISAalllA.;* 

12. 

ISYaI|ir*2 

13. 

is=imi ( ibA.isri 

14. 

JIHafl 

15. 

[50=15 

16. 

IE      (ISO-Hl      10.10.5 

17. 

3 

I5u=ISo-M 

IB. 

II 

I  1=11 

19. 

15 

CALL    CAi.C 

20. 

N/=i| 

21. 

NSHZ 

22. 

IE     lU»Jf-l)     25,2s,2u 

23. 

20 

A5=A(im-|  ) 

2*. 

rs  =  Y (n-i  I 

25. 

AlN-1 )=' IN) 

4b. 

r  (n-i  )»y  (im) 

27. 

A<N| B<S 

2B. 

r  In) =rs 

29. 

25 

ISalS.l 

30. 

Jl  =  l  1 

31. 

J'i 

IE     (  Is-'V)     4D.35. 35 

32. 

35 

lb=[S-4 

33. 

4  0 

||,A»INA  (  15) 

34. 

ll'»=l  *1  1  IS) 

35. 

iA^a [A. 1 OA 

3b. 

1 YialY. lur 

37. 

J|H=JlH»i 

38. 

IE     (JTH»|5(ill)     5li.-io.4S 

39. 

4s 

*hI  IE     l^lhl    Cv 

40. 

rff  MMN 

22 

?4 
2*i 

2« 

30 
3? 
3* 
3  6 
3i 


c 

50 

c 

5? 

c 

54 

C 

5* 

c 

SH 

c 

60 

c 

6? 

c 

64 

c 

66 

c 

64 

c 

70 

c 

1? 

c 

74 

r 

76 

c 

7H 

c 

HO 

c 

82 

c 

H4 

c 

46 

c 

tH 

c 

HIM 

c 

92 

c 

<V4 

r 

96 

c 

9H 

c 

1  1)0 

c 

1  02 

c 

1  OA 

c 

1  1)6 

c 

1  f)M 

c 

1  10 

c 

1  \? 

c 

1 1* 

c 

I  16 

c 

1  1" 

c 

1  20 

c 

122 

c 

1  24 

c 

1  ?6 

c 

1  2H 

c 

1  30 

r 

1  i? 

c 

1  34 

c 

1  36 

c 

1  3"t 

c 

1  4  ' 

c 

142 

r 

1  44 

c 

1  46 

c 

14* 

c 

1  so 

c 

152 

c 

1  S4 

r 

ISS 
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*9, 

50  • 

51  . 
52. 
S3. 
5*. 
55. 


9]  . 
^  . 
SI3. 


Si 

If 

(  ISSI       'St  /S 

SS 

3S 

If 

I  I  <-  1  Al.  I     SIS 

hi)  iV) 

bK 

If 

I  1  1-  1  Hi  1     SJS 

6S.yi 

Ion 
11.5 

lin 

1  h 


C 
12P 


L*S 

I  1" 


U'l 
US 
ISn 


ISS" 

c 


If     I  IS- J  SO)     'vSt  7u  i'Vt 
call    f.OLC 

so   r>>  isu 


if 
if 
if 
if 

if 
if 

IICP* 

IK     ( 

If  ( 
l«=l 
R'A  = 

f  r=? 

CALL 
I  A=  I 

do    r 


i    bdil*iliHO 

-Ml)      HS.HS.l*!) 

I     9ritl40tS)0 

-k.1)     S>S.SISil<.0 

AM(H?tIY2)  )     10o»  10I1.1S 

».^.I(;t«»2-ll  Jo5tl?5tlu5 

U»iUI««»ll»?iI>l*»M(I*ilt?)»AHlI<^iI»ill/».0 

-cv)    isiiin.no 

(1S-11)     U3.12u.115 


,»5 


I  r=  I  r ♦  i  ir 
IDYs-Illi' 

f  r  =  /.ii 
OLL    CALL 

i  t  =  i  r ♦  i ,j r 

If      (4M(M-],lY)-CV)      Ul).  US,1  JS 

...f  =   If  ♦  1 

*f  C  I  if)  =100»IA»IY 

IS=IS.4 

|  A  =  I  A? 

1  (  =  [  IV 

(,()    fn    3.1 

a  i  =  i  r 

if      (AM(iA-liIY)-CV)     14S.lSfi.lSM 

Nf  =     l*-»«  1 

ftFC  I   if)  =10(I»IA*IY 

l.Ll  -hS=  .Ll'.tb'l 
1)0     1SS    *=1,N 

A  («.)  =A  I  <.  I  -f  INN 

y (K) =Y(K)-KMZ 

CUNflNUt 


C    Sin, 

C 


t  C'wn   ui> 


I  Aft  ) 


U     l.j)     |/i/ 


» I7ui  lbi) 

•JtCL 

IA(1)  i  r  (  I)  t  la  1  till 

HI)  I     Cu 

-s)     (a(  1 1  .r  1 1 ) . 1  =  1  |N) 
i  .EU.ii)    ij.j    in    1 7 n 

I  .11  I  ."IAAA)  »»a*i|i  ( 1 1 
)  .(il  iMAA  Y  I  MAAiraY  (  I  I 
I  .LI  .Ml.JA)     MIN««A  (  I  ) 

r .minv)    mi <r  =  r  1 1 1 
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r)      LISTING    OF    SUBROUTINE    CALC 


Sooi  ijuTltJt    C»LC 


|C|     ,J,   „f    1HM      Iulr.rit'lJLllT  IOi   S     "[TmIn     daTa 

Cm  ,  IO'ih    L  I  nl     CO'J-vdI  'A  1Kb. 


■  I  A     b'iv.nt  ^      I"     LOC..  It 


IMF    oCAHb    StSTfcM 
kUwOut    OiMlWfcMSirr 


A.     Kfc  I  In    Iu-imKm, 
CUlL     tNuIr.FtnlMii 

ji  iiE     IS»n    . 


IS..<-t     — 

CALL    Cai.C 


•SCHIPflun    Of     COimThOL    CnKDS    - 

IH.jK     HtOUlHEO    -ALL     HtCtSbAKY     kAHA^fcTFHi    f»bSH)     VIA     uLAIvK     CO"M()l 


'JU  Ii<K  I)     1U11MIUI  I  ut  ■ 

-    -Jl)   it 


i)k  I  [U  j«l    ■iUHMIUI  Iiit  b    — 

-     IjUtit 


DIHFVIMOU    AMdOUf  1QU)  •     HFClHOO).     A  (  1  S  fl  0  )  <     rll'lUO).     lkl(J,3>.      1NM" 
lli     MY(h)i     lUH(2uli    *L(M),     YLilu),    /Llini,    i.Ab(lOiO),    QA(2),    ore*' 

^1.      '1/(2)1      H  I    (  1  .J  I   .      CIHKIAI4I,      S|H£|AU)«      n|!>l|«)l       1111*1.      .UKMIUJI, 

II 
coi<i.-(on   -.t  I'll  in 1 1 1*,  iy,  iua,  tor ,  i Si 1 1 1 ,  i</ , "(P,k,  jr  ,Pt  ,«£  c.cv  ,  a,  r .  ikr  , 

I  IfvA.  I  JV  .l)L  .AM,  IOP.Uk.UO.  ILnfliXL  i  YL.liL  iL".1.0»  ,ur  ,0/  i  *f"A  <  .Ml  ,SCAl.t,Y 
><MAX,ClHt  l«,5lMETAiFi'i2iFMNil;NIM,CL»l3iNLlivirS,HC:uHVi/,  I  1 1  ,  «UKK 

I  l  =  n 

u=ij.i    ■ 

IK     (I|)X»«2«II)Y»»2-1)     2(l»bi20 

U     I  ll)X)     m,  lu.  lb 

«  (  N  I  a  I  x 

/-lx 

[Y2=IY» |UY 

ljr=II)Y 

r(N)*((.iMlIXllY)-CV)/lAM(I»iIY)«AM(IXiIYH)|)»UY»/ 

HI-  I  OWN 

r  (i.i  =i  Y 
.=  1  < 

I1A=|DX 

[*^=  [  H»  IIjA 

(IN)i((a'lUAil()-tVI/(JH(I>ilYI»«»(Uililrll)«U*t« 
rd-  I'J'In 

I  A^=  I  A.  (||A 

1  r^s  I  r ♦  i  ut 

«=IX 

^  =  IY 

IJA  =  I')X 

■  J  Y=  1 1>  Y 

IHV»<  AMI  [A,  IY>  «A<-U  U2«  IY)  fAMlIx  1 1  Y2)  »ArM  lAH.  IY2)  > /4.0 

U      (  P  Y->»  .  II I     2S.jl).?b 

If     IDO-CV)     3u.Jii.J0 

«L='M(  Ia,1  ri-UCk 

i/  =  .S»  I  AL»Urk-CV)  /  ^L 

A  (Hi  =V«   U», 

(lui=««  >r ./ 
pr=  i  .ii 

HI-    I0><  .J 

»L  =  A  1  I  I  <2l  I  Y2I  -llCH 

v  =  .S»  (Ai.  •UCP-CV  I  /»L 

A  I  '-I  I  =  -  1/  » I  I A  ...  U  A 

r  (ui  =-i/»ur .  /«i)y 

HH    I  OH  J 

(■nil 


11 

100 

f) 

10? 

n 

1  04 

n 

1  0* 

n 

1  0* 

n 

1  10 

0 

)  1? 

n 

1 14 

n 

1 1* 

D 

1  1H 

o 

l?n 

0 

12? 

n 

12* 

n 

l?* 

n 

I?H 

n 

1  )<1 

n 

1  t? 

n 

1  34 

IJ 

1  Jft 

0 

1  3H 

n 

1 4') 

n 

14? 

Q 

1  s  - 

n 

\   I,f, 

i 

I  4  -i 

n 

1  41 

o 

IS? 

n 

1  S4 

ii 

1  ■j'i 

C-83 


o 

/l 

UJ 

y 

z 

'J 

§ 

i 

o 

X 

£ 

3 

. 

cr 

(\ 

o 

u. 

£>  **  f-    t>  i^  r-   jo 

CO 

(V  (\i  (\J 

o 

z 

H 

co 

_) 

UJ 

o 

9 

z 

o 

UJ 

~* 

LT 

" 

UJ 

^ 

r-  cr   4    -   ri  .t  it- 

U. 

-    rv  i\    <\    r» 

UJ 

^ 

cr 

CO 

T 

CO 

u 

T 

O 

-1 

I 

or 

u 

; 

o 

'.' 

7 

uf    -   cfi        u    •     a. 
—  —  (V  rv  f   n 

H 

— 

t 

< 

*.' 

K 

' 

^ 

CO 

~) 

I/) 

cr 

1/) 

O 

1/1        i/  y 


I         l  r 


O        'I  I  I 


Ti        l         i         r 


C-84 


•  ,.i  c 


CV 


1' 
or 
\)i 
Km? 
FiMh 

Hi 
IDA 
I  LIT 

ILhH 
[HA 

[NY 
I  OP 

I' 

lib 
I  I 

I  J 
!  « 
l*<- 


I  r, 
jr 


Ml 

vLII-.fcS 


SCM  E 
blnf  la 

I  I  I 


It  K  1 1 

1  IIU 

>  , 

?• 

en 

>• 

c 

J  . 

CD 

i. 

C 

1. 

c 

3. 

c 

). 

c 

3. 

c 

i. 

CI) 

?. 

Ci.i 

?. 

CO 

). 

c 

1. 

c 

s  , 

co 

la 

c 

). 

c 

). 

c 

<". 

CI) 

?. 

Cli 

?• 

Co 

,?. 

CI) 

)• 

C 

f. 

c 

1. 

c 

1. 

c 

c 

CO 

c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c< 
c 

CD 
CD 


CD 

CI) 

CU 

C 

CD 

CD 

C 

CD 

CI 


3. 
12, 


lt>. 

V. 

12. 


■itl-cl-E    CI-  b 


Jo. 


1 1. 

JO. 
I*. 

lo. 


3  I . 


in.  ?».  31 . 

12.  S5.  iS, 


3'.  38.  3S>. 

U*  3o. 

11.  31.  3H. 

I*,  it.  34. 

12.  ii.  32. 


if 

STCrful  is 

nt-F  i  >t U 

4 

MM  ll    Uf 

NUMrlr*     (' 

f     kEF  F  nfc  ,Ct  S 

l/l 

IOIAI 

NUM4>H    (1 

F    «Ki  E*t  »rt  s 

r1  •  1 

vh    '       it- 

I         ,-H 

y     »tl- 1  <h    'f.F  •■ 

C-85 
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